Singular perturbation techniques in the gravitational self-force problem 
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Much of the progress in the gravitational self-force problem has involved the use of singular per- 
turbation techniques. Yet the formalism underlying these techniques is not widely known. I remedy 
this situation by explicating the foundations and geometrical structure of singular perturbation the- 
ory in general relativity. Within that context, I sketch precise formulations of the methods used 
in the self- force problem: dual expansions (including matched asymptotic expansions), for which I 
identify precise matching conditions, one of which is a weak condition arising only when multiple 
coordinate systems are used; multiscale expansions, for which I provide a covariant formulation; 
and a self-consistent expansion with a fixed worldline, for which I provide a precise statement of 
the exact problem and its approximation. I then present a detailed analysis of matched asymptotic 
expansions as they have been utilized in calculating the self-force. Typically, the method has relied 
on a weak matching condition, which I show cannot determine a unique equation of motion. I 
formulate a refined condition that is sufficient to determine such an equation. However, I conclude 
that the method yields significantly weaker results than do alternative methods. 

PACS numbers: 04.20. Cv, 04.25. Nx 
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I. INTRODUCTION 

Perturbation theory is a venerable field of study in gen- 
eral relativity (GR). In fact, because of the complexity 
of the Einstein field equation (EFE), most physically rel- 
evant analytical results in GR rely on perturbing away 
from a known solution. However, most of the founda- 
tional work in this area has focused only on descriptions 
of regular perturbation problems — problems in which a 
regular power-series expansion yields a uniform asymp- 
totic approximation to a true solution. The underlying 
formalism of regular perturbation theory in GR has been 
studied extensively [1-3], and it has been shown that 
any regular asymptotic expansion of the field equations 
yields a pcrturbative solution that approximates at least 
one exact solution, at least locally [4-6]. 

However, many physically interesting systems must be 
modeled as singular perturbation problems — problems in 
which a regular power-series expansion fails to provide a 
uniform asymptotic approximation. Indeed, one of the 
most successful areas of research in GR, post-Newtonian 
theory, centers on a singular perturbation problem. As 
one would expect, the foundations of that particular 
problem have been studied extensively [7, 8], and it is 
now known that there exist a large class of exact solu- 
tions possessing post-Newtonian expansions [9-12]. But 
general discussions of singular perturbation theory in GR 
arc lacking; Kates has provided the only such discus- 
sion [13], and his emphasis was on providing an overview 
of the geometrical structure of singular problems, forego- 
ing any discussion of particular methods. 

Along with the post-Newtonian expansion, another 
problem of great significance is also singular: the mo- 
tion of an asymptotically small body through an exter- 
nal spacetime. This problem will be the primary focus 
of the present paper. Study of the point particle limit 
is less well developed than that of the Newtonian limit. 



and it has typically focused on proving that at leading 
order, a small, uncharged body behaves as a test mass, 
moving on a geodesic of the external background space- 
time (see, e.g., Refs. [14-19]). However, in recent years, 
the advent of gravitational wave astronomy has enjoined 
a need to go beyond the test particle approximation. 
Specifically, more accurate approximations are required 
to model extreme mass-ratio inspirals (EMRIs) , in which 
a compact body (such as a neutron star or black hole) 
of mass m ~ Mq spirals into a supermassive Kerr black 
hole of mass M ^ (1O'*-1O^)M0 lying at the center of a 
galaxy; see Refs. [20, 21] for an overview of these systems. 
Extreme mass-ratio inspirals are a potentially important 
source of wave-signals for the planned gravitational- wave 
detector LISA [22], with a predicted rate of several to 
one hundred detectable events per year [23, 24]. 

For an EMRI, an expansion in the point particle limit 
roughly corresponds to an expansion in powers of the 
mass ratio m/M ^ e. (For the remainder of this pa- 
per, I assume all variables have been scaled by a global 
lengthscale such as M, such that I can write, e.g., m ^ e, 
where e is dimensionless.) At leading order in a regular 
expansion, the small body behaves as a point particle and 
moves on a geodesic of the spacetime of the large body. 
At sub-leading order, the metric perturbation generated 
by the body exerts a force on it, the dissipative part of 
which causes it to very slowly spiral into the large body. 
The acceleration of this inspiraling worldline, caused by 
the body's interaction with its own gravitational field, is 
called the gravitational self-force. A general, analytical 
expression for the self-force in an arbitrary vacuum back- 
ground spacetime was first derived by Mino, Sasaki, and 
Tanaka [25] and Quinn and Wald [26] ; that expression is 
now known as the MiSaTaQuWa equation. 

This is a singular perturbation problem for two rea- 
sons. First, near the body, the first-order metric per- 
turbation is that of a point particle, behaving as a delta 
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function with support on the particle's worldline. If the 
body is anything other than a black hole, the error in 
the approximation (i.e., the difference between the exact 
and approximate metrics) will then be unbounded in a 
small neighborhood of the body; and even in the case of a 
black hole, the expansion will fail, since the second-order 
Einstein tensor will contain products of delta functions 
and hence be ill-defined as a distribution. Therefore, the 
approximation breaks down due to rapid changes, on the 
lengthscale ^ e, near the body. The second reason the 
problem is singular is the existence of secular errors (i.e., 
ones that accumulate over time) . If we use the first-order 
solution with a point particle source, then the Bianchi 
identity constrains the particle to move on a geodesic; 
at higher order in the expansion, this geodesic worldline 
is corrected by small deviation vectors [27, 28]. But the 
true path of the small body spirals into the large black 
hole, eventually deviating by a very large amount from 
the leading-order, geodesic worldline; generically, the er- 
ror in the regular expansion will be unbounded on any 
time-interval [0, p > 0. Therefore, the approxima- 
tion also breaks down due to slow, cumulative changes 
over the radiation-reaction timescale trr = (In an 

EMRI, this is the time required for the particle's energy 
and angular momentum to undergo an order- 1 change, 
since their rate of change is proportional to the self-force, 
itself of order e.) 

In singular perturbation theory, in order to overcome 
these types of errors, one makes use of general expan- 
sions,^ of the form 

JV 

/(a;,£) = ^e"/(")(x,£) + 0(£^+i). (1) 

n=0 

Unlike a regular power-series expansion, here the coeffi- 
cients f^'^\x,e) are allowed to depend on e; but they are 
nevertheless constrained to be of order 1, in the sense 
that there exist positive constants k and Eq such that 
|/(")(a;,e)| < fc for < e < Eq, but lime^o /^"^(a;, e) ^ 
(unless /'"'(a;,£) is identically zero). Put simply, the 
goal of a general expansion is to expand only part of a 
function's e-depcndence, while holding fixed some spe- 
cific £-dependence that captures one or more of the func- 
tion's essential features. There are two common types 
of general expansions: composite expansions and multi- 
scale expansions. The first of these is designed to over- 
come the failure of a regular expansion near some sub- 
manifold (such as the position of the small body in an 
EMRI) , while the second is designed to overcome secular 
errors (such as the deviation of the true motion from the 
leading-order, geodesic approximation). 

Composite expansions are patched together from a fi- 
nite number of regular expansions. For example, in the 
EMRI problem, if we use a rescaled radial coordinate 



^ In Rcf. [28], I instead used the term "singular expansions". 



f = r je near the body, then a regular expansion at fixed 
f could be accurate on the scale r ~ er, where an ex- 
pansion at fixed r fails; this inner expansion can then 
be combined with an outer expansion valid for r ~ 1, 
yielding a composite expansion that is uniformly accu- 
rate both near and far from the body. I will refer to a 
pair of inner and outer expansions such as this as dual 
expansions. One can make use of dual expansions in a va- 
riety of ways, but historically, they have been used most 
often in the method of matched asymptotic expansions, 
in which the two expansions are first partially determined 
in their respective domains of validity, and then any re- 
maining freedom in them is removed by insisting that 
they agree in a domain of mutual validity.^ They were 
first used in fluid dynamics to analyze the behavior of 
a low-viscosity fluid near a boundary. In the context of 
GR, since the pioneering work of Burke [29], who stud- 
ied the effect of radiation-reaction on a post-Newtonian 
binary, and D'Eatli [16, 17], who studied the motion of 
black holes, dual expansions have mostly been utilized 
for two purposes: determining waveforms by matching 
wave-zone expansions to near-zone expansions (see, e.g., 
the review [30]), and determining equations of motion 
for small bodies by matching an inner expansion near 
a body to an outer expansion in a larger region (see, 
e.g., [15, 31-36]). The latter method, in various forms, 
has often been used to derive the MiSaTaQuWa equation 
[25, 27, 28, 37, 38]. (See Ref. [28] for a more thorough 
review.) 

Composite expansions are suitable only when the dif- 
ferent lengthscales dominate in different regions — e.g., 
the metric varies on the short lengthscale ^ e near the 
small body, while it varies on the background lengthscale 
~ everywhere else. In contrast, multiscale expansions 
are suitable in situations where multiple lengthscales are 
relevant everywhere in the region of interest. These ex- 
pansions begin directly with the generalized form (1), 
where the coefficients /^"■\x,e) depend on some specific 
function of e — for example, in a two-timescale expansion, 
/("^ is assumed to depend on a time coordinate t and a 
slow-time coordinate t = et, which allows the expansion 
to capture both short-term and long-term effects. Re- 
cently. Hinderer and Flanagan [39] have suggested a two- 
timescale expansion of the EFE and equation of motion 
in the EMRI problem; the method has also been applied 
to more restricted expansions of self-forced equations of 
motion [40, 41]. 

In addition to these standard expansions, there has 



^ My nomenclature is not standard. In most of the literature on 
singular perturbation theory, the term "matched asymptotic ex- 
pansions" has the meaning that I have here assigned to dual 
expansions. However, I have opted to follow the usage in re- 
cent literature on the self-force, which has used the term to refer 
specifically to the procedure in which the field equations in the 
inner and outer expansions are first solved and then the solutions 
are made to match. 



3 



been one further general expansion utilized in the self- 
force problem: a self-consistent expansion with a fixed 
worldline. In this expansion, rather than allowing a 
functional dependence on an e-dependcnt function from 
spacetime to R, such as a slow-time coordinate, one al- 
lows a functional dependence on a function from M to 
spacetime. This allows one to consider a metric per- 
turbation that is a functional of an exact, e-dependcnt 
worldline, bypassing the constraint that the worldline 
must be a geodesic at leading order. This idea underlies 
the methods of post-Newtonian theory, and it has been 
assumed in some form in much of the literature on the 
self- force, including in the earliest derivations of the MiS- 
aTaQuWa equation [25, 26]. In Ref. [28], this approach 
was first formalized in terms of a systematic approxi- 
mation scheme. That scheme makes use of two general 
expansions: an inner expansion accurate near the small 
body, and an outer expansion accurate in the external 
background spacetime. In the outer expansion, the met- 
ric is treated as a functional of the exact worldline 7, 
and it is then expanded while holding that dependence 
fixed: g = g + eh^^^[-/] + e^h^^lij] + .... See Ref. [28] 
for further details. Like a two-timcscale expansion, this 
fixed-worldline approach promises to eliminate the sec- 
ular errors of a regular expansion; however, the precise 
relationship between the two methods remains to be ex- 
plored. 

In this paper, I seek to accomplish two goals. My first 
goal is to extend the geometrical description of singular 
perturbation techniques. The geometrical picture of reg- 
ular perturbation theory in GR is well known; the anal- 
ogous description of singular perturbation theory, pro- 
vided by Kates [13], is less well known. In Sec. Ill, I 
review Kates' work and then go beyond it by describ- 
ing in some detail the three techniques used in the self- 
force problem: dual expansions, multiscale expansions, 
and the self-consistent expansion presented in Ref. [28]. 
I focus on making precise statements pertaining to these 
methods, with the aim of clarifying previous work. With 
regard to dual expansions, I identify two matching condi- 
tions that may be used in the method of matched asymp- 
totic expansions — a strong matching condition analogous 
to the one used in traditional singular perturbation the- 
ory, and a weak matching condition that arises when 
working in multiple coordinate systems. With regard to 
multiscale expansions, I provide a covariant description. 
With regard to the self-consistent expansion, I extend 
the discussion of Ref. [28] by providing an exact formu- 
lation of the self-force problem and a brief discussion of 
the gauge freedom in the expansion. 

My second goal is to provide a precise formulation 
of the derivations of the gravitational self-force using 
matched asymptotic expansions. Section IV consists of a 
new version of this derivation, along with detailed analy- 
sis and discussion. Matching was used in some of the ear- 
liest derivations of the MiSaTaQuWa equation [25, 37]. 
Because it can determine the first-order equation of mo- 
tion from a first-order outer expansion, it is technically 



far less involved than the method of Rcfs. [15, 27, 28], 
which determine the nth-order equation of motion by 
solving the field equations to n + 1st order in a region 
around the body. However, as I will discuss, the method 
as it has been utilized in the self-force problem has re- 
lied on the weak matching condition, which is too weak 
to actually determine an equation of motion. In order 
to determine an equation of motion, additional assump- 
tions must be made. I pinpoint these assumptions and 
formulate a refined matching condition. Because of the 
variety of assumptions required, I conclude that matching 
yields weaker results than one would expect from tradi- 
tional singular perturbation theory; and it yields weaker 
results than those obtainable through the more laborious 
approach of Refs. [15, 27, 28]. The method can surely be 
improved, the number of assumptions reduced, but such 
improvement may prove more difficult than using alter- 
native methods. 

Before proceeding to these two goals, I begin in Sec. II 
with a review of traditional singular perturbation meth- 
ods in applied mathematics. Such a review is warranted 
both to develop the concepts required in the later sections 
and to make the basics of the theory more widely known. 
More detailed overviews of the subject can be found in 
numerous textbooks (e.g., Refs. [42-46]). Among these, 
the text by Kevorkian and Cole [44] covers the broadest 
range of topics, and the text by Eckhaus [46] provides 
the most rigorous treatment. 

II. TRADITIONAL SINGULAR 
PERTURBATION THEORY 

I begin by defining some useful notation. First, I define 
the following order symbols: for x € K", 

• f{x, e) = 0{({e)) if there exist positive constants k 
and e* such that |/(a;,e)| < /c|C(£)l for fixed x and 
< e < e*. 

• f(x, e) =^ o(C(e)) if lim = Q at fixed x. 

./(x,e) = 0,(C(£)) if f{x,e) = 0{C{e)) and 
/(x,e)^o(C(e)). 

For example, 5e + 2e3/2 = 0(e) = 2e^^'^, 5e + 2e3/2 ^ 
o(e) = 2^3/2, and be + 2e^/^ = ^(e) ^ 2£3/2. 

In general, we are concerned with the asymptotic be- 
havior of functions, rather than the behavior of functions 
evaluated at particular locations. This means we need a 
norm appropriate for a function. Also, a central issue 
in perturbation theory is whether or not an approxima- 
tion is uniformly accurate in a region of interest, where 
uniformity is defined as follows: 

• f{x,e) = 0{(^{e)) uniformly in a region D if 
there exist positive constants k and e* such that 
\\f{x,e)\\D < k\C{e)\ for < e < e*, where 
II ■ \\d = sup^g^ I • |. 
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Analogous definitions hold for o and Og- These defini- 
tions provide a more useful measure of the asymptotic 
behavior of a function. 

Finally, I define several important asymptotic quanti- 
ties: relative to an asymptotic sequence {Cni^)}, where 
C„-|-i(£) = o(C„(£)), 

• f{x,e) is an A^th-order asymptotic approximation 
of fix, e) = 0,(1) if fix, e) - f (.t, e) = o(CAr(e)),3 

• fix,e) is an iVth-order asymptotic solution to a 
differential equation D[f(a::,e)] = if D[/(x,£)] = 
oiCNie)), 

N 

• fix,e) = Y,Q^ie)f^^\x,e), where f^^\x,e) = 

0^(1), is an A^th-order asymptotic series. \i f'^"-^ is 
independent of e, then the series is said to be reg- 
ular (sometimes called Poincare-type) ; if not, then 
it is general. If, in addition, / is an asymptotic 
approximation to f , then it is an A^th-order asymp- 
totic expansion of f . 

The most common asymptotic sequence is {£"}, which I 
will use almost exclusively in this paper. Note that we are 
typically uninterested in whether or not an asymptotic 
series converges as A^ -> oo. In fact, even if / is an 
asymptotic series that both converges and asymptotically 
approximates f, the function that it converges to might 
not be f . 

In any given perturbation calculation, one almost al- 
ways calculates an asymptotic solution to an equation. 
Determining whether or not an asymptotic solution is 
also an asymptotic approximation to an exact solution 
is typically far more difficult. It is, however an essential 
step in proving the reliability of an expansion, because 
an asymptotic solution may not be an asymptotic ap- 
proximation to an exact solution (more pathologically, 
an asymptotic approximation may not be an asymptotic 
solution [42]). 

General asymptotic expansions are a powerful tool for 
solving singular perturbation problems, which are defined 
by the failure of a regular expansion to provide a uniform 
approximation. This failure is often signaled by a change 
of character in the governing differential equation when 
e — 0: for example, a hyperbolic equation might degen- 
erate into a parabolic equation. The inaccuracy of a reg- 
ular expansion is also frequently signaled by its failure to 
satisfy a given boundary condition, or by the expansion 
growing without bound in a system that we have reason 
to believe should be bounded. Typically, the underly- 
ing origin of the failure is the presence of distinct length 
scales, one of which appears only for e > 0. Often, this 
means that the exact solution to a problem is singular 



^ In the case f (a;, e) = Os(ffc(e)), functions would be rescalcd by 
Cfe(e) before making comparisons. 



at e = 0. Hence, in singular perturbation problems, we 
assume that e G (0,e*], which allows us to take the limit 
e ^ 0, but which prevents us from setting £ = 0. In this 
section, I consider two types of systems: first, systems in 
which the exact solution undergoes a rapid change near a 
submanifold; second, systems in which rapid changes oc- 
cur throughout the region of interest. In the first type of 
system, dual expansions can be used to construct a uni- 
form general expansion; in the second type, a multiscale 
expansion can be used. 

Before proceeding, I define two final pieces of notation. 
-0* and tp* denote, respectively, the push-forward and 
pull-back corresponding to a map ip. So, for example, if 
/ is a function of coordinates x, and x = ?/'(x), then ip^f 
is the function rewritten in terms of x. <I>^ f denotes the 
Ath-ordcr regular asymptotic expansion of f in the limit 
of small £, holding fixed the coordinates of which f is a 
function. 



A. Dual expansions 

Dual expansions are typically used to solve bound- 
ary value problems in which the solution exhibits rapid 
change in a very small region (or a finite number of such 
regions). The regions of rapid change are referred to as 
boundary layers. Frequently, this rapid change prevents a 
regular expansion from satisfying a given boundary con- 
dition, though a "boundary layer" can sometimes arise 
away from any boundary. The usual means of solving 
these problems is to make use of two regular expansions: 
an inner expansion /;„ that is expected to be valid in 
the boundary layer, and an outer expansion /out that 
is expected to be valid outside of it. Suppose we have a 
one-dimensional problem with coordinate r, and that the 
boundary layer is at r = r?, and has a thickness ^ ^ie). 
Then the outer expansion is simply a regular series at 
fixed r, and the inner expansion is a regular series at 
fixed values of the rescalcd coordinate 7~ = (r — rf,)/i5(£); 
this can be written as : r i-^ f. The inner expan- 
sion allows us to capture changes over the Icngthscale 
Sis), since r is of order unity when the original coordi- 
nate r is of order Sie). Note that if we treat the prob- 
lem on a two-dimensional plane with coordinates ir,e), 
then the inner and outer expansions can be visualized as 
expansions along flow lines deflncd by r ^constant and 
r/e =constant, as shown in Fig. 1. 

For simplicity, suppose that ri, ~ and that boundary 
data is given at r = and r ~ I. In this case, the 
outer expansion typically fails to satisfy the boundary 
condition at r = 0, but it can be made to satisfy the 
condition at r = 1; conversely, the inner expansion can 
satisfy only the condition at r = 0. This leaves each of 
the expansions underdetermined. The basic idea of dual 
expansions is to fully determine them by insisting that 
they agree in some region of mutual validity. Suppose 
that /in is an A^th-order asymptotic approximation of 
f in a region Din, and /out is an Ath-ordcr asymptotic 
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FIG. 1: Upper: inner limit (dashed curves) and outer limit 
(dotted curves) in the (r, e)-plane. Lower: the same limits in 
the (f, e)-plane. The inner limit is defined by e — >■ 0, f = r/e 
fixed; the outer limit, by e — ^ 0, r fixed. From the perspective 
of the inner limit, the outer limit sends all points to infinity 
if — !■ oo). From the perspective of the the outer limit, the 
inner limit sends all points to zero {r Q). 



approximation in a region Dout • Then by definition, f — 
/in = o{Cn{£)) in An and f - /out = o{Qn{£)) in Aut- 
Subtracting the second equation from the first, wc have 
the overlap matching condition: 



/out - /in = o(CAr(e)) in Dout n Z^i, 



(2) 



Note that this condition relies on the existence of the 
overlap region Dout H Din- If that region is empty, then 
the condition is vacuous. And it may not be obvious that 
such a region ever exists, since the inner expansion triv- 
ially appears to be valid only for f ^ 1, and the outer 
expansion only for r ~ 1. However, if /out('') is a uniform 
asymptotic approximation to f (r) on an interval [a, 1] for 
constant a, then it is also a uniform approximation on 
the extended interval [Cz(e), 1] for some (i{e) = o(l); sim- 
ilarly, if /in(^) is a uniform approximation to f (f) on [0, 6], 
then it is a uniform approximation on the extended in- 
terval [0, 1/Cj(e)] for some Ci(e) = o(l) [45, 46]. If we 
have access to the exact function f, then we can explic- 
itly determine the overlap of these extended regions. But 
in a typical application, without access to an exact solu- 
tion, we must make use of the overlap hypothesis, which 
states that the overlap region exists. In order to imple- 
ment the overlap matching condition, one then simply 
assumes that the constructed asymptotic series of a given 
order are asymptotic approximations of the same order. 



and one then takes the overlap region to be the region in 
which /out - /in = o(CAr(e))- 

In this paper, I will not make direct use of the over- 
lap matching condition. Instead, I will use a second, 
simpler matching condition, which I will refer to as the 
coefficient- matching condition: 



(3) 



In this matching condition, we match results term-by- 
term in the expansions. On the left-hand side we have the 
inner expansion (<I>™?/;*f) written as a function of r (via 
"0*) and then expanded in the outer limit; on the right- 
hand side, we have the outer expansion (•I'^f) written as 
a function of f (via -0,) and expanded in the inner limit, 
and then rewritten as a function of r and re-expanded. 
(The right-hand side requires an extra expansion in or- 
der to remove terms that would appear as higher-order 
terms in the inner expansion. Refer to Appendix A for 
an illustrative example.) We can write this schematically 



(4) 



meaning that when the inner expansion is re-expanded 
for small e at fixed r, and the outer expansion is re- 
expanded for small r at fixed e, the two results must 
agree term by term. We can then, for example, equate 
coefficients of e"r™ on the left- and right-hand sides. 

If we define the buffer region by the inequalities e ^ 
r <|; 1,** this equation states that the inner and outer 
expansions must agree term by term when they are both 
expanded in the buffer region. In other words, if the exact 
solution is expanded first for small e at fixed r/e (yielding 
an inner expansion), and then expanded at fixed r (or in 
other words, for r ^ e), it must agree, term by term, 
with the result of expanding first for small e at fixed r 
and then expanding for r <^1. 

From the perspective of the inner limit, the buffer re- 
gion lies at asymptotic infinity (f — > oo); from the per- 
spective of the outer expansion, it lies asymptotically 
close to r = 0. From this we can intuit a still sim- 
pler matching condition, which I will call the asymptotic 
matching condition: 



where fl^'^ and /oJ^i are the leading-order terms in, re- 
spectively, the inner and outer expansions. 

The three matching conditions I have discussed are 
obviously related. In fact, one can derive the asymp- 
totic matching condition and (a condition similar to) the 



lim /j 



f(o) 



(0) 



hm /-^"^ 
inn /out, 



(5) 



* In the applications of matched asymptotic expansions in GR, the 
meanings that I have assigned to the terms "overlap region" and 
"buffer region" are often conflated, and the terms are often used 
interchangeably. For the sake of clarity, I distinguish between 
the two. 
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coefficient-matching condition from the overlap hypothe- 
sis. However, one should realize that the overlap hypoth- 
esis is merely sufficient to arrive at those two matching 
conditions: it is not necessary. Functions exist that do 
not satisfy the overlap hypothesis but nevertheless satisfy 
the coefficient-matching condition, for example [46] . 

Once a matching condition has been used to fully de- 
termine the inner and outer expansions, one can con- 
struct a composite expansion that is uniformly accurate 
on the full domain of the problem. This expansion con- 
sists of the sum of the inner and outer expansions, minus 
the terms that are common to both in the buffer region. 
Explicitly, 

/comp = + r'i'Ty^*^ - ^^rv^^'&rf, (6) 

which we can write schematically as 

/comp /out ^ ./in ^r/out- ('^) 

Note that this is a general asymptotic expansion of the 
form E„Cn(e)F(")(r,r/5(£)). 

The matching conditions presented here can be used 
in a variety of ways. Most traditionally, in the method of 
matched asymptotic expansions, the inner and outer so- 
lutions are determined as far as possible using the govern- 
ing differential equation and boundary conditions, and 
then they are fully determined by imposing a match- 
ing condition. However, one can instead begin with the 
matching conditions to restrict the general form of one or 
the other of the expansions, or to develop a general ex- 
pansion in the buffer region. In Appendix A, I illustrate 
the method of matched asymptotic expansions with a 
simple example. For further information on dual expan- 
sions, see Refs. [45, 46]. 

B. Multiscale expansions 

In some systems, rather than a rapid change occur- 
ring near a submanifold, rapid changes occur over the 
entire region of interest (in other words, the "boundary 
layers" are dense in the region). When studying these 
systems, one cannot form a uniform approximation by 
making use of regular expansions in separate regions and 
then combining them. Instead, one must assume a gen- 
eral expansion from the start. 

Suppose for simplicity that the rapid changes occur 
on a scale ~ 1 and the slow changes occur on a scale 
~ and that we seek an approximation to f(t,e) 

that is uniform on the time- interval [0,1/e]. Then we 
proceed by introducing a fast time variable (j) = </'(^: £) 
and a slow time variable i = i{t, e) satisfying ^ = uj{t, e) 



^ One could rescale the variables such that the rapid changes occur 
on the scale ~ e and the slow changes on the scale ~ 1 , to accord 
with the description of a region dense with boundary layers. 



and ^ = euj{t,e), where ui and Cj are uniformly Os{l); 
changes in are of the same order as changes in t, while 
t changes appreciably only after t changes by a very large 
amount. (In the simplest case, we have (j) = t and t = 
et.) We invert the transformation in order to write the 
frequencies as functions of the slow time alone: w(t, e) 
and a)(i, e). I next note that while setting a; = 1 will lead 
to large errors on a timescalc 1/s (consider, for example, 
attempting to approximate cos(i + et) by cost), setting 
CJ = 1 will lead to large errors only on extremely long 
timescales outside our range of interest. Hence, I will 
take the slow time to be given by t = et. The remaining 
frequency, w, must be determined over the course of the 
calculation. To make such a goal feasible, I assume lo 
possesses a regular expansion X]n>o 

I next assume that f (t, e) can be written as a function 
F[(f), i, e), and that F possesses a regular expansion: that 
is, 

f{t,e) - F((/),i» = ^Cn(£)^^^"^(0,t). (8) 

n 

Suppose f is to satisfy some differential equation D[f] = 
0. After making the substitution f{t,e) = F{(f),i,e), we 
use the chain rule to convert derivatives with respect to 
t into the sum of partial derivatives ^ = Lo{i,e)-^ + 

e-^. We then arrive at a partial differential equation in 

terms of (f> and t. Now, the essential step in a multiscale 
expansion consists of treating (j) and t as independent 
variables at this point; that is, -F is taken to be a solution 
to the PDF for arbitrary values of 4> and t. Given this 
assumption, after substituting the expansions for F and 
Lu, we can solve the equation by setting the coefficient of 
each C„ to zero. If I did not assume that F solves the 
equation for arbitrary (f> and t, then the e-dependence 
scattered throughout (j){t, e) and i{t, e) would prevent us 
from solving the equation in this manner. 

Treating and t as independent is equivalent to work- 
ing on an enlarged manifold with coordinates {(f>,t,e). 
The solution manifold on which f lives is a submanifold 
defined hy cp = <p{t,e) and t = i{t,e). (See Fig. 2 for an 
illustration in the simple case where 4> = t and t = et.) 
Determining w can be viewed as a step in determining 
this submanifold; in fact, we can note that the transfor- 
mation from the extrinsic coordinates a;" = (^, e) to 
the intrinsic coordinates = (^j s) defines a set of basis 
vectors e" on the submanifold, given by ef = (u;,e, 0) 
and = (8,(1), t, 1). 

Because we are provided with sufficient boundary data 
for an ODF, rather than for a PDF, we must place 
some constraints on the function F. The most com- 
monly imposed constraint is the non-secularity condition, 
which says that integration constants must be chosen 
such that any secularly growing term vanishes. Other 
possible constraints include the demand that each co- 
efficient F^''\4>,i) is a periodic function of and the 
demand that each coefficient be unique. Obviously, one 
must apply such constraints judiciously and systemati- 
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FIG. 2: In a multiscale expansion, we work on a manifold of 
larger dimension than that on which the original problem is 
posed. The solution is eventually evaluated on the submani- 
fold defined hy t = et, shown here in grey. 

cally. 

Of course, this procedure relies on a host of assump- 
tions. There is no guarantee that the exact solution f 
possesses an asymptotic expansion of the form (8); and 
even if it does, there is no guarantee that the terms in the 
expansion wih necessarily solve the equation for arbitrary 
(j) and t. However, this method is extremely successful 
in practice. If one instead assumes a regular expansion, 
then the dependence of f on t will be expanded in pow- 
ers of et. These powers of t will eventually grow large, 
such that terms initially supposed to be high order be- 
come as large as the lower-order terms, preventing the 
expansion from providing a uniform approximation. In 
many cases, one can avoid this secular growth by using 
a rigorous method of averaging, which removes the rapid 
time dependence and recovers only the leading-order slow 
time-dependence. However, if one requires the fast time 
dependence as well, then the two-time method offers the 
most powerful means of doing so. 

In Appendix B, I illustrate the above ideas with an 
example. In that example, I demonstrate that even if 
the assumptions of the multiscale expansion fail at some 
order, the lower order solution can still yield a uniform 
asymptotic approximation; in addition, the failure of the 
assumptions is made manifest over the course of the cal- 
culation. For further information on multiscale expan- 
sions, see Ref. [44]. 



C. Singular versus regular perturbation theory 

We should now take note of the essential differences be- 
tween regular perturbation techniques and singular per- 
turbation techniques. When a regular expansion of an 
exact solution f is substituted into a differential equation 
D[f] = 0, the coefficients in the expansion are guaranteed 



to solve a hierarchy of differential equations, simply by 
setting the coefficients of each power of e to zero. Hence, 
when constructing a regular series solution to a differen- 
tial equation, one can determine each term in the solution 
solely from the given differential equation (and its atten- 
dant boundary conditions). But a general expansion of 
an exact solution is not guaranteed to satisfy any such hi- 
erarchy, because the coefficients in the expansions depend 
on e. Hence, when constructing a general series solution 
to a differential equation, one must impose some extra 
conditions upon it — e.g., satisfying the overlap hypothe- 
sis in the method of matched asymptotic expansions, or 
satisfying a PDE rather than an ODE in the method of 
multiple scales — which are not guaranteed to be satisfied 
given only the form of the general expansion. 

This means that proving general properties of solutions 
is much more difficult using singular perturbation theory. 
In regular perturbation theory, one can construct proofs 
of the form "given an exact solution to such and such a 
boundary value problem, if it possesses a regular asymp- 
totic expansion then that expansion has such and such 
behavior" ; in singular perturbation theory, we must ap- 
pend further hypotheses to this statement. However, if 
we seek a very strong statement about the solution to 
a problem, we must in any case go beyond the form of 
such a proof: we must also prove that the exact solu- 
tion actually does possess the assumed expansion. This 
is a difficult feat regardless of whether the assumed ex- 
pansion is regular or general. While it is usually easier 
in the case of regular expansions, techniques do exist for 
handling singular perturbation problems (see Ref. [46] 
for examples). Furthermore, general expansions provide 
asymptotic solutions where regular series cannot, and 
they provide uniform asymptotic solutions. Hence, in 
most cases of interest, their advantages far outweigh any 
disadvantages. 



III. PERTURBATION THEORY IN GENERAL 
RELATIVITY 

In GR we typically do not begin with a predetermined 
manifold with predetermined boundary conditions that 
uniquely determine an exact solution. Instead, the man- 
ifold is (mostly) determined by the leading order "back- 
ground" solution to the Einstein equations. Within that 
manifold, we define boundary conditions that uniquely 
determine the perturbations. This somewhat complicates 
the problem, but it also makes the assumptions of sin- 
gular perturbation theory more reasonable: since we do 
not seek an approximation to a unique exact solution to 
a given boundary value problem, but only an approxima- 
tion to some exact solution to the EFE, it is eminently 
reasonable to impose the supplementary conditions re- 
quired to construct general expansions. 
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A. Regular perturbation theory 

Before describing singular perturbation theory in GR, 
I will briefly review regular perturbation theory. In its 
most geometric description, the formalism begins with 
a 5D manifold M carrying a 5D metric g"^" of signature 
(0, — , +, +, +) and a non-negative scalar field £ : A/" ~> M. 
The manifold is foliated by 4D submanifolds A^g defined 
by de = 0, such that M ^ Me x R. When restricted to 
act on dual vectors tangent to M.^, g"^" can be inverted, 
inducing a 4D Lorentzian metric ge. Each member of the 
family of metrics {ge} is taken to be an exact solution 
of Einstein's equation at fixed e. A regular expansion of 
the pair (A^e, ge) is an expansion around a known "base" 
pair (A^Ojff = go)- This expansion is performed by first 
defining a one-to-one relationship between points on A^g 
and A^e via a diffeomorphism (^e : A4o — > A/Je, called an 
identification map. This map induces a flow on M with a 
tangent vector field X that is non-vanishing and nowhere 
tangent to the submanifolds Me- (Note that one could 
instead begin with a vector field and derive from it an 
identification map, but beginning with the identification 
map will be more useful in formalizing general expan- 
sions.) 

In this context, the regular expansion ge = .g + 
^^^i is given by an expansion along the flow in- 

duced by X : 

0:(ge)=e^^-g|^^, (9) 

where 4>*ege is the pull-back of ge onto the base manifold 
A/Io, and £x is the Lie derivative along the vector X. 
The background metric and the perturbations of it have 
clear geometrical interpretations: the background metric 
g is the restriction of ge to the submanifold defined by 
£ = 0; the first-order perturbation eft,'^^' = £(»Cxg)|^_Q is 
the product of the "distance" e along a flow line and the 
rate of change of g in the direction of the flow; and so on. 

A choice of gauge corresponds to a choice of idcnti- 
flcation map ipe- A different choice, say ipe, leads to a 
different tangent vector held Y , which in turn leads to a 
change 

reige) - reige) = {e^^ " e^^'-)g|_^^. (10) 

By expanding the exponentials, one flnds that this in- 
duces changes /i^"^ — > /i'"^ -I- A/i'"^. At first and second 
order, the changes are given explicitly by 

A/id) = £^,,,5, (11) 

A/.(2) = i(i;e,,,+4^^)g + i:c,„/iW, (12) 

where ^(i) = {Y - X)\^^^ and $(2) = [A:,r]|^^g are vec- 
tor fields in the tangent bundle of M.q. Note that ^(1) 
and ^(2) are linearly independent, so they can be chosen 
independently. In terms of coordinates, they correspond 
to the near-identity transformation 

^ ^ _ ^ + 1,2 (efi),,efi) - ^(-2)) + 0{e^) , 

(13) 



where the components on the right hand side are in the 
original coordinates cc" . We say that the vectors are 
the generators of the gauge transformation. (See Ref. [3] 
for the precise meaning of this phrase.) 

B. Singular perturbation theory 

Although singular perturbation techniques have been 
utilized in many calculations in GR, the only formal de- 
scription of them was provided by Kates [13]. I will re- 
view his description in this section, before extending it in 
the following sections. A singular perturbation problem 
is characterized by the limit e — > being singular: ge may 
not exist at e = 0, the topology or dimension oi M.e may 
change between e = and e > 0, etc. This means that 
the 5D manifold N does not in general contain a "base" 
manifold A4o; instead, it is given by A/" ~ Me x (0,e*]. 
Hence, one cannot generically build an approximation on 
the limiting manifold by finding derivatives of the exact 
metric at £ = 0. Instead, one works on a "model mani- 
fold" Mm, on which one constructs a family of approxi- 
mate solutions 

gM{e)^g{e) + Y,e^^h(-\e). (14) 

n>l 

The topology of the model manifold is taken to be com- 
patible with the leading-order metric g{e). If there exists 
an identification map ipe : Mm — > Me, which maps a 
region Um C Mm to a region Ue C A^ej such that ffM(e) 
uniformly approximates (flge in the region Um as £ — > 0, 
then (?m(£) is a uniform asymptotic approximation (as 
measured in some suitable norm) to the exact solution 
in the region Ue ■ Once again, the identification map in- 
duces a family of curves in the 5D manifold Af, but these 
curves will not, in general, continue smoothly to a base 
manifold Mq. 

As an example, consider a post-Newtonian expansion, 
which is singular [7, 36, 47, 48]. The Newtonian limit 
is given by£ = i;/c— >0; where v is the supremum of 
the velocities in the system. In this limit, the light cones 
of the spacetime fold out into spatial surfaces, and the 
time-components of the metric blow up — alternatively, 
if we consider the inverse metric, we see that its time 
components vanish, such that it degenerates into a 3D 
spatial metric. Hence, we can infer that the manifold de- 
fined by £ = corresponds to the 3D spatial manifold of 
Newtonian theory.® In this case, the model manifold and 
background metric are taken to be those of Minkowski 
spacetime. 

In the next two subsections, I formulate dual expan- 
sions and multiscale expansions within this framework. 



The singular nature of the Newtonian limit is also signaled by 
the fact that hyperbolic wave equations become elliptic Poisson 
equations as the speed of gravity's propagation becomes infinite. 
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The description of dual expansions follows that given by 
Kates [13], which built on the work of D'Eath [16, 17]; 
however, I more carefully formulate the matching con- 
ditions, specifically stressing the distinction between the 
strong matching condition (used in traditional singular 
perturbation theory) and a weak matching condition (of- 
ten used in GR). My discussion of multiscale expansions 
is original to this work. 

In the final subsection, I formulate the self-force prob- 
lem as a free-boundary value problem. I then discuss 
means of solving the problem within the context of a 
self-consistent expansion with a fixed worldline. 



C. Dual expansions 

Though most of the description in this section carries 
over to a more general situation, I will restrict it to the 
pertinent case of a family of exact solutions containing 
a body of mass ^ e, on a family of manifolds M.^. In this 
section, I only sketch the formalism of inner and outer 
limits for this system; in Sees. IIIE-IV, I provide further 
discussion of concrete applications of these limits. 

Suppose that we are given two coordinate systems on 
M.^: a local coordinate system X" = {T,R,Q^) that 
is centered (in some approximate sense) on the small 
body, and a global coordinate system x". For exam- 
ple, in an EMRI, the global coordinates might be the 
Boyer-Lindquist coordinates of the supermassive Kerr 
black hole (though we could consider the case in which 
both coordinate systems are centered on the small body) ; 
the local coordinates might be Schwarzschild-type coor- 
dinates for the small body. The local coordinates cover 
some region Dj around (and possibly inside) the body, 
while the global coordinates cover a larger region De out- 
side the body. Assume, without loss of generality, that 
the two coordinate systems have overlapping domains, 
and that they are related by a map (j)^ : .t" i— >■ in the 
region D = Dj Ci De- 

A regular outer expansion gE{x,e) = g{x) + eh'-'^'> (x) + 
...is constructed by taking the limit £ —J- at fixed x". In 
this limit, the body shrinks toward zero size as all other 
distances remain roughly constant. For simplicity, I as- 
sume that this limit continues to a base manifold Mq. 
However, the limit certainly does not exist on a remnant 
curve 7''°-' corresponding to the "position" of the small 
body — for example, if one takes a regular limit of the 
Schwarzschild metric in Schwarzschild coordinates, then 
there is no limit defined at coordinate values correspond- 
ing to r = 0. Hence, I take the model manifold in the 
outer expansion to be TW^ = Mq U 7^"-', and I take the 
external background metric to be g = go. Of course, this 
construction is not essential; the model manifold need not 
be defined by setting e = in this way. But at the very 
least, for the outer expansion to be regular, we require 
g = lim^^o g. 

A regular inner expansion gi{X,e) = gB{X) + 
eH^^'>{X) + ... is constructed by taking the limit e — ^ 




FIG. 3: Regular inner limit (dashed curves) and outer limit 
(dotted curves) on the manifold Af ~ A4s x (0, e*]. The inner 
limit is generated by the map ipj from the manifold Mi on 
which the interior background metric of the isolated small 
body lives; these curves terminate at a point {T — To,R — 0) 
on 7'*^^ The outer hmit is generated by the map (pE from 
the manifold Me on which the external background metric 
lives. The external manifold in this case is taken to be equal 
to Mo U 7^'^'. The gray region is a surface of constant R, 
which converges to the £-independent worldline 7'"-' . 



at fixed values of the scaled coordinates X" = ip{X°') = 
((T - To)/e,i?/e,e^). This limit is naturally singular: 
it follows flow lines that converge at a single point de- 
fined by (T = To,i? = 0) in A^^;. ExpUcitly, since 
the metric written in these coordinates has the form 
g ^ e^gBapdX^dX^ , all distances vanish at £ = 0. As 
discussed by D'Eath [17] (see also Ref. [27]), to make the 
limit regular, one must use the conformally rescaled met- 
ric ge = £~^ge. This rescaling effectively "blows up" the 
distances in spacetime, such that as £ — > 0, the size of 
the small body remains constant while all other distances 
arc sent to infinity; thus, the inner limit serves to "zoom 
in" on a small region around the body. The background 
spacetime defined by £ = is then defined by the metric 
gs of the isolated small body, and the approximation is 
built on a model manifold Aij with the topology of that 
spacetime.^ 

The outer and inner expansions are related to the exact 
solution via identification maps ipE : Me ^ Me and 
ifi : Mi —> Me, which respectively identify points on 
Me and Mi with points on Me- (See Fig. 3.) These 
two identification maps induce a map (p : Me ^ Mi, 
given hy cj) = ipJ^oipE, which has the identical coordinate 
description as the original transformation between the 
global and local coordinates. Gauge transformations in 



^ Note that Mi gcncrically differs from Me- Consider, for exam- 
ple, the case of a small black hole orbiting a large black hole. 
The manifold Mi possesses a singularity at the "position" of 
the small black hole but is otherwise smooth, while the manifold 
Me possesses a singularity at the "position" of the large black 
hole but possesses a smooth worldline where the small black hole 
should be. 
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the outer and inner expansions are generated by vector 
fields <^"(a;) and ^"(X), which take their respective values 
in the tangent bundles oi Ais and Aii- Note that a 
gauge transformation in the outer expansion generically 
corresponds to a finite coordinate transformation in the 
inner expansion, due to the rescaling of the coordinates. 

A uniform composite expansion is formed on a model 
manifold Ai ^ Me by cutting out a portion oi A4e and 
stitching part of Mi into the excised region. The local 
and global coordinates each cover a patch of M, identi- 
fied with the patches Dj and De on Me via the maps 
(pi and (pE- The uniform metric is constructed on this 
manifold by adding together the inner and outer approx- 
imations in each coordinate system, then removing any 
"double-counted" terms that appear in both metrics. 

How would one go about constructing such a uniform 
approximation if one did not have access to an exact so- 
lution? Just as in traditional perturbation theory, one 
would construct two separate asymptotic solutions to 
Einstein's equation, but now in possibly two different co- 
ordinate systems and on possibly two different manifolds. 
If one assumes that the two asymptotic solutions arc ap- 
proximations of a single exact solution, and if there exists 
an overlap region on the manifold M in which both ap- 
proximations are valid to the same order, then they must 
agree in that overlap region. In this case, "agreement" 
is defined by the existence of the unique map (p that re- 
lates the two expansions. As usual, we will not worry 
about a specific overlap region, but instead expand the 
two solutions in the buffer region. 

However, before performing that expansion in the 
buffer region, one must write gi and g^; in the same coor- 
dinate system. Let us choose this system to be the local 
coordinates X". Then, adapting Eq. (3), the matching 
condition reads 

(15) 

On the left, we begin with the exact metric in the local 
coordinates A"". It is then expanded to ?Tith-ordcr in an 
inner expansion, by transforming into scaled coordinates 
X" via t/i (along with an appropriate conformal rescaling) 
and expanding. Next, it is expanded in the buffer region 
by re-expressing it in the unsealed local coordinates and 
expanding to fcth order; this is equivalent to an expan- 
sion of the inner solution for e. On the right-hand 
side of the equation, we begin with the exact metric in 
the global coordinates . It is expanded to mth order in 
those coordinates, yielding an outer expansion. It is then 
transformed to the starting point of the left-hand side, 
by transforming to the local coordinates via 0, then to 
the scaled local coordinates via ip^ then re-expanding to 
?7ith order to yield an inner expansion. Finally, it is ex- 
panded in the buffer region by transforming back to the 
unsealed local coordinates and re-expanding. The con- 
tent of this equation is that the expansion in the buffer 
region must be the same whether it is obtained by first 
performing an outer expansion or by first performing an 



inner expansion. Schematically, we can write 

$e5/(^) = *i?0*.9£, (16) 

which states that if the inner and outer expansions are 
written in the same coordinate system, then they must 
yield the same expansion in the buffer region. 

But this is decidedly not the matching condition that 
has been used in practice. Instead, what has been done 
in practice is the reverse: first, expand the two solu- 
tions in the buffer region, and only afterward find the 
coordinate transformation between them. This is accom- 
plished by setting up a second local coordinate system 
= 0-y(x") = {t,r,6^) centered on a worldline 7 in 
Me', for example, these might be Fermi normal coordi- 
nates, and in the case of regular expansions, they would 
be centered on 7^°'. The outer expansion is then written 
in these local coordinates and expanded for small r, un- 
der the presumption that r ~ i?. After performing this 
expansion (and the expansion of gi in the buffer region), 
one seeks a unique transformation 0buf : X°' that 

maps the buffer-region expansion of gE into the buffer- 
region expansion of g/.^ Schematically, we can write 

<^e9l{X) = 0buf*$r07*5i5- (17) 

On the left, the inner expansion gi is expanded in the 
buffer region in the local unsealed coordinates X". On 
the right, the outer expansion g^; is transformed to the 
local coordinates Y" via <j>^, then it is expanded in the 
buffer region (i.e., for small r). Hence, the two buffer- 
region expansions are written in two different coordinate 
systems: the inner expansion in the coordinates X", and 
the outer expansion in the coordinates Y". So, in order 
to make a comparison, as the final step on the right-hand 
side, the buffer-region expansion of gE is transformed to 
the coordinates X" via 4>huf*- In short, Eq. (17) states 
that if gi and gE are expanded in the buffer region, then 
the resulting expansions must be related by a coordinate 
transformation. 

I will call Eq. (15) the strong matching condition and 
Eq. (17) the weak matching condition. The weak condi- 
tion follows from the strong condition, but not vice versa, 
and one can easily imagine situations in which the weak 
condition would be satisfied while the strong condition 
would not. In the weak matching condition, because 
the metric is already expanded for small r before (f>huf 
is determined, 4>hui will itself be written as an expan- 
sion. Thus, the weak matching condition only requires 
an asymptotic approximation of 4>hui (or, equivalently, of 
= (phuf^P'y)- Of course, one can only ever determine an 
asymptotic approximation — but in the strong matching 
condition, the approximation is for small e, rather than 
for both small e and small r. This essentially reduces 



One can see that if everything is correct, the various transforma- 
tions must be related as </> = (/)buf ° 07 • 
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(phui to a gauge transformation in the buffer-region ex- 
pansion defined by R (or r) and e both being smaU. As 
mentioned above, a gauge transformation in the outer ex- 
pansion corresponds to a finite coordinate transformation 
in the inner expansion, and vice versa. Hence, any choice 
of gauge on Me must be compatible with the choice of 
background coordinates on Aii (and vice versa). The 
two matching conditions insist on this compatibility to 
differing extents. 

One should note that though the description in this 
section makes use of two regular expansions, as in tra- 
ditional matched asymptotic expansions, the same gen- 
eral description holds for two general expansions. The 
only differences arc that g ^ lim^^o g and that there is 
no need to conformally scale the metric to arrive at the 
inner expansion. This is particularly important if one 
wishes to allow the internal metric to vary on its "nat- 
ural" timescale T ~ 1 (i.e., the timescale determined by 
the mass of the small object). If the metric near the 
body varies on this timescale, then in the unsealed co- 
ordinate time T = eT, the metric will have a functional 
dependence on the combination T/e, which will be singu- 
lar in the limit e — > 0. Thus, if both the expansions are 
to be regular, the internal metric can vary only on the 
external time T, corresponding to an internal slow evo- 
lution depending only on eT. In other words, regularity 
requires that the internal solution varies quasistatically 
(see D'Eath's discussion [17]). Of course, for e > 0, one 
could construct a general inner expansion that is iden- 
tical to the regular inner expansion by rescaling R only, 
instead of both T and R, and then simply assuming the 
inner expansion varies quasistatically; using this method, 
a global-in-time expansion can be constructed, and the 
metric is never conformally rescaled. Also, by using this 
method; one can remove the quasistatic assumption en- 
tirely. 

Finally, before moving to the next singular perturba- 
tion technique, I will note that just as in traditional 
singular perturbation theory, there is a distinction be- 
tween what I have called the overlap region and the 
buffer region. The buffer region corresponds simply to 
e ^ R 1. In order for us to express the outer solution 
in terms of the field R, the buffer region must lie within 
the region D, where both the local and global coordinate 
systems apply, but the size of the region is independent of 
the order of accuracy of the inner and outer solutions. As 
discussed in Refs. [15, 27, 28, 33, 36], one can extract con- 
siderable information about the metric — and in particu- 
lar, equations of motion for the small body — by working 
entirely within the buffer region, without ever construct- 
ing explicit inner and outer solutions or making use of 
an overlap hypothesis. This information is typically ex- 
tracted by defining the mass and current moments of the 
body in the buffer region, which is possible because the 
buffer region lies at asymptotic infinity from the perspec- 
tive of the inner expansion. Solving the Einstein equation 
then determines the evolution of these moments. In par- 
ticular, an evolution equation for the body's mass dipolc 



informs us of the motion of the body's center of mass 
relative to the chosen local coordinate system, providing 
an equation of motion for the body. 



D. Multiscale expansions 

In the method of multiple scales, changes on both 
short and fast time scales occur throughout the space- 
time. Thus, one cannot construct a uniform asymptotic 
approximation based on combining only two limit pro- 
cesses. If we consider a two-timescale expansion, with a 
fast time t and a slow time t = et, there are only two 
limits that can be easily envisioned: the slow-time limit 
e — >■ at fixed t, which follows a congruence of curves 
in J\f that tend toward t oo as e — ?> 0; or the fast- 
time limit e at fixed t, which follows a congruence of 
curves that tend toward t — > 0. However, in a multiscale 
expansion, both quantities are to be kept fixed. As dis- 
cussed above, this is accomplished by treating them as 
independent variables. 

Consider the case of an expansion that holds fixed both 
a set of coordinates and some scalar field ({x, e) sat- 
isfying ^ - o(l): 

g(x,£)=g(x,C) + Es"'^^"H^,C)- (18) 

n>l 

In the simplest case, C is equal to the product of e and 
one of the coordinates. When substituting this multi- 
scale expansion into Einstein's equation, one would treat 
C and as independent coordinates on an extended, 5D 
manifold Aie] these 5D manifolds are stacked atop one 
another to form a 6D manifold A/" A^e x M^. The limit 
e — is taken at fixed values of both C, and x, and the 
actual solution is obtained by restricting the expansion 
to the submanifold defined by C = C{x,e). 

As in traditional perturbation theory, one might re- 
quire a fast-time variable cj) that differs from the given 
coordinate time. Indeed, one might use any coordinates 
one likes on Aie- To provide some flavor of the expan- 
sions, in this section I will define gauge transformations 
and sketch a multiscale expansion of the EFE for the 
simple case with coordinates (a;",C); I will occasionally 
provide details given the additional simplifying assump- 
tion 5pC = £V^{x, C) for some = Os{l). 

Note that the gauge transformations in this expansion 
differ from those of a regular expansion. Gauge transfor- 
mations arc generated by transformations of the form 

x°'^x'°' =x°'~eC{x,0+O{e^), (19) 

where f = Os(l)- In order to determine the effect of this 
transformation, I expand the Lie derivative as 

£^ = £P+£^^\ (20) 
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where, e.g., for a vector ^(x, C) and a tensor T''^(x, C), where Rl^pvcr{xX) is constructed from g and V^''^ 



(21) 



dC dxP ^ " dC dxP ' ^ ^ dC dx" ' 

(22) 



These definitions are independent of the behavior of C,. 
But in the particular case that = ^^(^JjOi the 
gauge transformation generated by a vector e^(a;, C) can 
be written as 

A/i(i) = (23) 
A/i(2) ^ i£W^W^^ ^ j^W^ii) + ^a)^^, (24) 

Similarly, I expand the covariant derivative as 

\V%x, C) = (V^") + V^i)) V^ix, C) (25) 

where is compatible with g{x,C,{x,e)), ^^f^ is com- 
patible with g at fixed ^, and X^''^' is compatible with g 
at fixed x. Explicitly, 



Vi°V''(x,C) 



I no : 



(26) 



vj^^^''(-'C) = ^^ + rW;^,y'', (27) 

where the Christoffel symbols are given by 



(28) 



t-i(1)q: 1 

2 5 \ g^j + g^fS g^ g^S 



(29) 



The "correction" X^^^'' ensures that the total covariant 
derivative \^ is compatible with the C-depcndence of g. 
Note that all three derivatives are metric compatible: 
\7g = 0, V(°)g = 0, and V^^^g = 0. 

By writing the Lie derivative in terms of the covari- 
ant derivative, we can express the gauge transformation 
generated by a vector field e^(x, C) as 

= 2V((")^,) (30) 

+ rV,("'/.L^j + Vj°)r + 2V(«e« (31) 

assuming that — eV^{x, C). 

Note that because the background metric g depends 
on C, the Riemann tensor constructed from it can be 
expanded in powers of e: 



Rfipvcrix-, C) — ^upucri-^^ C) + ^-^lipiycr (^^j C) 



ppva 



(32) 



R\ipvcr{x, C) contains one V*^^-' derivative, and Rppua{x, C) 
contains two V*-^-* derivatives. The nth-order perturba- 
tion of the Ricci tensor can be similarly expanded as 

5''Rp^[h] = EL=o£"'^"-Rp"^[^]- This means that the 
vacuum Einstein equation Rpi, — becomes 

R^p} = 0, (33) 
<54°)[/.W]=4i), (34) 
<54°)[/.(2)] ^ _ _ 5^R^)[hi% (35) 



Similarly, the Bianchi identity on the background. 



: 0, becomes 
V(o) . G^") = 

y(0) . ^ . 

^p ^Up ^p ^Up 5 

^p ^Wp ^p ^Up 5 



(36) 
(37) 
(38) 
(39) 



where a dot indicates contraction over and v. And the 
Bianchi identitity on the full spacetime, g^'^^VpG^p = 0, 
can be expanded schematically as 

V(o) .5GW[/i(i'] =0 (40) 



,5r(")[/i(i)]. (g(i) + 5gW[Mi) 

V(o) . j2^(o)[/i(i)] (41) 



where SV'-"^'^ [h^^^] 



l„a5/y(0), (1) 



'^s^^'^ ^^p-y) ^ ^^'^ the leading-order Einstein equation Rp'J = 
and the Bianchi identity gP'^VpGyp = have already 
been imposed for compactness. 

Note that the C-depcndence of g allows the background 
to slowly react to the perturbation. Determining this re- 
action is the backreaction problem, which has been stud- 
ied extensively in the past. 

Recently, Hinderer and Flanagan [39] have constructed 
a significantly more complicated two-timescale expansion 
tailored to EMRIs. In their method, all dynamical vari- 
ables (i.e., the metric and the phase space variables of 
the worldline) are submitted to two-timescale expansions; 
this expansion captures both the fast dynamics of orbital 
motion and the slow dynamics of the particle's inspiral 
and the gravitational backreaction on the background 
spacetime. Since the metric and the worldline are re- 
lated by the EFE, it is assumed that the metric can be 
written as a function of the phase space variables of the 
worldline. On each timeslice, the limit e — >■ is then 
taken with the phase space variables held fixed. Specif- 
ically, the true worldline is specified by a set of action- 
angle variables {J{t,e),Lp{t,e)) and a slow time variable 



(0) 
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t. Expanding for e — )■ with tp and t held fixed resuhs 
in a sequence of fast-time and slow-time equations. In 
the fast-time equations, wherein t (and therefore J) is 
treated as a constant, the metric is a function of ip only; 
in other words, it is a functional of the geodesic that is 
instantaneously tangential to the true worldline. From 
this it follows that the leading-order fast-time equation 
yields a metric perturbation sourced by that geodesic, 
as in regular perturbation theory. However, that is only 
at fixed t — the true worldline and metric perturbation 
emerge by allowing the variables to vary with t, with a 
t-dependence determined from the slow-time equations. 

Although exceedingly useful for EMRIs, this proce- 
dure relies on the background metric being stationary 
at fixed t, such that it has no fast time dependence, and 
on the geodesic motion in that background being inte- 
grable, such that the metric can be written in terms of 
the action-angle variables. In addition, it requires one to 
determine the slow evolution of the background metric, 
which has not yet been done. 



E. Self-consistent expansion 

In Hindcrcr and Flanagan's formalism, the metric is 
written as a function of the phase space variables on the 
worldline, and then both the metric and those variables 
are submitted to a two-timescale expansion. The formal- 
ism I will now describe is a generalization of this: the 
metric is written as a functional of the worldline, and 
then the metric is expanded with that worldline held 
fixed. In order to motivate this approach, I will first 
provide a formulation of the exact problem to which wc 
seek an approximate solution. 

We wish to determine the mean motion of a small, spa- 
tially bounded matter distribution. (For the moment, I 
neglect the case of a black hole.) In principle, we have 
some matter field equations to go along with the EFE 
for this blob of matter. As governed by the field equa- 
tions, the boundary of the blob traces out some surface 
in spacctimc. In the interior of the boundary, the matter 
density is finite, and in the exterior it vanishes. To de- 
termine the motion of the body, we seek the equation for 
the generators of this boundary. This is a free-boundary 
value problem [49], in which some boundary values are 
specified on a boundary that is free to move. In the con- 
text of bodies in GR, this problem has received some 
study [50, 51], but it is still far from understood, and it 
must certainly be tackled numerically. 

To make progress with an approximation scheme, I re- 
formulate the problem. I surround the body by a tube F 
embedded in the buffer region, such that for e ^ 0, the 
radius of the tube vanishes. For the moment, consider F 
to be defined by constant radius R = TZ{e) in the local 
coordinates X". I assume that the body is fairly widely 
separated from all other matter sources, such that outside 
of F there is a large vacuum region fl. I also assume that 
F is in vacuum; since it lies in the buffer region around 



the small body, this means that I must restrict my ap- 
proximation to a small body that is sufficiently compact 
to not fill the entire buffer region. Now, since the tube 
is close to the small body (relative to all external length 
scales), the metric on the tube is primarily determined 
by the small body's structure. In other words, the infor- 
mation about the body has now been transplanted into 
boundary conditions on the tube. Recall that the buffer 
region corresponds to ^ — > oo. Hence, on the tube, we 
can construct a multipole expansion of the body's field, 
with the form . I assume that the local coordi- 

nates X" are mass-centered, such that the mass dipole 
term in this expansion vanishes. (See Ref. [52] and ref- 
erences therein for discussion of multipole expansions in 
GR; see Refs. [33, 52] for discussion of mass-centered co- 
ordinates in the buffer region; see, e.g., Ref. [53] for fur- 
ther discussion of definitions of center of mass.) This, 
then, is another free-boundary value problem: we must 
determine the equations of motion of the generators of 
the tube, given the boundary values of the metric on it, 
and in particular, given that the body lies at the "cen- 
ter" of it. With this formulation, we can also determine 
the motion of a black hole, rather than just a matter 
distribution. 

Now suppose that I want to represent the motion of 
the body through the external spacetime (g, Me)-, rather 
than through the exact spacetime. As we can see from 
Fig. 3, this is easily accomplished by using the regular 
limit and taking the motion to be represented by the 
remnant worldline 7^*^^. However, as mentioned in Sec. I, 
on long timescales this will provide a very poor repre- 
sentation of the motion. (See Ref. [28] for a detailed 
discussion.) 

Let us consider this from another direction. Assume 
that we were given the exact solution on A^g, along 
with the coordinate transformation 0^ between the lo- 
cal coordinates X" and the global coordinates in the 
buffer region. At fixed i? = TZ{e) = o(l), we could 
write this transformation as x" — (j)~^{T,TZ,Q^). In 
the limit of small e, TZ becomes small as well, mean- 
ing that this transformation can be expanded as a; = 
<f>-^{T,0,&^) +0(1), where 9^ is an arbitrary choice 
of angles. This transformation thus defines a curve 
z"(r,£) EE (j)j^{T,0,Q^) in the external manifold Me- 
Since the small body is centered "at" R = 0, this curve 
defines a meaningful long-term representative worldline 
7. If we expand (j)e{T, 0, 8^) for small e, then it will not 
provide a uniform transformation between the inner and 
outer coordinates; it will contain secularly growing errors 
of the form et. So, instead, in order to construct a uni- 
form asymptotic solution, when constructing the exter- 
nal approximation, one must hold 7 fixed. Determining 
7 then amounts to determining the "location" at which 
the (mass-centered) inner expansion is to be performed. 

Since we will never be seeking </) directly, and in case 
the inverse of (j)^ does not exist at i? = 0, allow me to 
present the final reformulation of the problem. Define 
a tube F£;[7] C Me such that it is a surface of con- 
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FIG. 4: Fixed-worldline expansion of a family of spacetimes. 
The dotted lines correspond to the outer limit, which lets 
the body shrink to zero size but keeps its motion fixed. The 
dashed lines correspond to the inner limit, which keeps the 
size of the body fixed. Here I display the singular inner 
limit, which does not rescale the inner time coordinate; hence, 
dashed lines originating at different times will terminate at 
different points on 7. The worldline lies in the manifold 
A4e = Mo but it does not correspond to the remnant 

curve 7'°-' defined by the regular limit; instead, it is allowed 
to have e-dependence, and it is determined by the particular 
value of e at which an approximate solution is sought. 

stant radius r in Fermi normal coordinates centered on 
a worldline j G A4e- Using the map ip e from the reg- 
ular expansion, this defines a tube T — (Pe^J'e)- Now, 
the problem is the following: what equation of motion 
must 7 satisfy in order for T to be mass-centered (in the 
sense that the mass dipole of the inner expansion vanishes 
when mapped to M.^ via ipi)l Note that the worldline is 
a curve in the external manifold AAe- It should not be 
thought of as a curve in the manifold AA^ on which the 
exact metric ge lives; in fact, if the small body is a black 
hole, then there is obviously no such curve. 

In order to determine the equation of motion of the 
worldline, I consider a family of metrics g£;(.T,£;7) 
parametrized by 7, such that when 7 is given by the 
correct equation of motion, we have 5£;(a;, e; 7(e)) = 
Lp*^g^{x). The metric in the outer limit is thus taken 
to be the general expansion 

g(a;, e) = gE{x, e; 7) = g{x) + h{x, e; 7), (42) 

where 

Ne 

hix,e-n) = 5^£"/i^")(x;7) +0(£^-+^) ■ (43) 

n=l 

Solving Einstein's equations will determine the worldline 
7 for which the inner expansion is mass-centered. I will 
call this either a self-consistent or a fixed-worldline ex- 
pansion. In it, the perturbations produced by the body 
are constructed about a fixed worldline determined by 
the particular value of e at which one seeks an approxi- 
mation. Refer to Figs. 3 and 4 for a graphical comparison 
between this expansion and a regular one. 



In the remainder of this section, I present a sequence 
of perturbation equations that arise in this expansion 
scheme, along with a complementary sequence for the in- 
ner expansion. The equations were originally presented 
in Ref. [28]. In that paper, I described a particular, 
lengthy method of solving the equations and deriving 
equations of motion. In Sec. IV, I will discuss an alterna- 
tive approach using the method of matched asymptotic 
expansions. 

My sequence of perturbation equations relics on a par- 
ticular choice of gauge. I discuss the gauge freedom in the 
self-consistent expansion, and the effect of gauge trans- 
formations on the equation of motion, in Appendix C. 

1. Field equations in outer expansion 

To begin, I surround the body with a worldtube T, 
where F is embedded in the buffer region, such that the 
field on it can be found from either the inner or outer 
expansion. I seek a solution in a vacuum region n outside 
of F; I further specify that il consists of the future domain 
of dependence of FUS, where S is a spacelikc initial-data 
surface. 

Now, recall that in a multiscale expansion, the ex- 
panded equations are solved by assuming that they are 
valid for arbitrary values of the slow-time variable t, not 
only on the true solution manifold defined hy t ~ et. 
Similarly, in the fixed-worldline expansion, one method 
of solving the expanded EFE will consist of assuming 
that it is valid for arbitrary worldlines; the true solution 
is found by choosing the true worldline. Solving the EFE 
with an arbitrary worldline seems to require reformulat- 
ing it in a "relaxed" form before expanding it, such that, 
for example, the linearized equation does not immedi- 
ately determine 7 to be a geodesic. To accomplish this, 
I assume that the Lorenz gauge can be imposed every- 
where in fl on the entirety of h, such that L^[h] = 0, 
where is the operator defined by 

L,[h] = - yX'Hf'p- (44) 

In Appendix C, I discuss the validity of this assumption. 

With this choice of gauge, the vacuum Einstein equa- 
tion R^i, = is reduced to a weakly nonlinear wave equa- 
tion that can be expanded and solved at fixed 7, leading 
to the sequence of wave equations 

E^4h^^^]=0, (45) 
E^4h^^^]=25^R^4h^^\ (46) 

where E^^^ is the relativistic wave operator defined by 

E^^Ah] = (^S'V'V, + 2i?/,'^) V- (47) 
More generally, we can write the Tith-order equation as 

£^..[4"^]=^^:'[4'^■•■>r'^7], (48) 
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where the source term S)Il consists of nonlinear terms 
in the expansion of the Ricci tensor. These equations 
can be solved for arbitrary T (and hence for abritrary 7). 
The formal solution is given by 

an 

j^G-P^'^'^Yl'dV. (49) 

where Gapa'P' is the retarded Green's function for E^p] 
I adopt the conventions of Ref. [54] for the Green's func- 
tion. This formal solution requires boundary data on V . 
Since T lies in the buffer region, the boundary data on it 
can be provided by the inner expansion, as discussed in 
Ref. [28]. 

In this scheme, the equation of motion is determined 
only once the perturbative EFE, rather than just the 
wave equation, is satisfied. One means of ensuring that 
the EFE is satisfied is to ensure that the gauge condition 
is satisfied. I assume that the acceleration of 7 possesses 
an expansion 

a^(i,£) = a(,o)(<)+£aW(i;7) + .... (50) 

This is an expansion of a particular function of time on 
the worldline; it does not suggest a multiplicity of world- 
lines with differing accelerations. Substituting this ex- 
pansion, along with that of h^^, into the exact gauge 
condition L^[h\ = and solving with arbitrary 7, we 
arrive at the sequence of equations 

L(°) [/.W] = 0, (51) 

L'^'[h^^]--LT[h'^l (52) 



where = L[/]L^„(o,, L^^^f] is linear in a^), 

L(2)[/] is linear in ' and quadratic in ' , and so on. 
More generally, for n > the equations read 

n 
m— 1 

In these expressions, L[f] is first calculated on an ar- 
bitrary worldline, and then to find L^'^\ the expan- 
sion of the acceleration is inserted — while still holding 
7 and fixed. These equations will determine suc- 
cessive terms ail^\ At the end of the calculation, when 
has been determined to the desired accuracy, Wfi] is 
evaluated for the particular worldline with acceleration 

0'^l — a{f''+£a^^'' + ..., just as at the end of a two-timcscalc 
expansion, the solution is evaluated for t ~ et. 

Beyond formalizing the self-consistent approach as a 
systematic approximation scheme, this method differs 
in one key respect from previous derivations of self- 
consistent equations of motion. Earlier derivations gen- 
erally result in equations containing derivatives of the ac- 
celeration; in the case of an uncharged body in vacuum. 



this is seen in the gravitational antidamping term discov- 
ered by Havas [55] (as corrected by Havas and Goldberg 
[56]). Such equations are unphysical, exhibiting, for ex- 
ample, runaway solutions. Traditionally, they have been 
made well behaved via an a posteriori "reduction of or- 
der" [26, 54]. However, my assumed expansion of the 
acceleration automatically yields a wcU-behavcd, order- 
reduced equation of motion. Furthermore, the expansion 
of the acceleration was necessary to split the gauge con- 
dition (or, equivalently, the Bianchi identity) into a se- 
quence of exactly solvable equations. Hence, we can see 
that a systematic approach, in which one requires exact 
solutions to the perturbation equations, eliminates the 
ill-behaved equations of motion that have plagued prior 
self-consistent derivations. 

In particular, this method differs from the gauge- 
relaxation procedure that has been used historically in 
the gravitational self- force problem [25-28] . In that pro- 
cedure, one constructs an approximate solution to the 
regular linearized Einstein equation by (i) solving the 
linear wave equation and (ii) ensuring the equation of 
motion enforces the relaxed gauge condition L^[eh^^'^] = 
O(e^). The relaxed gauge condition accomplishes the 
same goal as my reformulation of the Einstein equation 
into relaxed form: it allows the body to move on an ac- 
celerated worldline rather than the geodesic worldline en- 
forced by the regular linearized Einstein equation. How- 
ever, this procedure is an a posteriori corrective measure, 
rather than part of a systematic expansion, and unlike 
Eq. (51), it results in an ill-behaved equation of motion 
requiring further a posteriori correction in the form of 
order-reduction. 

Before proceeding to the inner expansion, I note that 
the first-order acceleration a*-^^ is determined by Eq. (52), 
which requires the second-order metric perturbation . 
Hence, if the gauge condition is to be used to determine 
the equation of motion, one must solve the second-order 
wave-equation in order to determine the first-order accel- 
eration. However, in the method of matched asymptotic 
expansions, the inner expansion at any given order gener- 
ically contains information that is of infinite order in the 
outer expansion. Therefore, if matching is used, it may 
be possible to derive an equation of motion using only 
the first-order term in the outer expansion. This will be 
the route explored in Sec. IV. 



2. Field equations in inner expansion 

For the inner expansion, I assume the existence of some 
local polar coordinates X°' = {T, R,Q^), such that the 
metric can be expanded for e — > while holding fixed 
R = R/e, Q^, and T. This leads to the ansatz 

g(x,£) = <7j(r,i?,e^,e) 

= gBiT,R,e^) + H{T,R,e^,e), (54) 
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where H at fixed (T, R, O"*) is a perturbation beginning 
at order e. The leading-order term gsiT, R, 6"*) at fixed 
T is the metric of the small body if it were isolated. For 
example, if the body is a small Schwarzschild black hole 
of ADM mass em(T), then in Schwarzschild coordinates 
gB{T,R,&'^) is given by 

ds^ = -(l- 2m{T)/R^ dT^ + " 2m{T)/R^ \^dR^ 

+ e^Ff {dQ^ + sin^ 6^$^) . (55) 

Since the metric becomes one-dimensional at e: = 0, the 
limit £ is singular. As mentioned in Sec. IIIC, the 
limit can be made regular by rescaling time as well, such 
that T = (T — Tq) / e, and then rescaling the entire metric 
by a conformal factor 1/e^. This is equivalent to using the 
above general expansion and assuming that the metric gs 
and its perturbations are quasistatic (evolving only on 
timescales ^ 1). Both are equivalent to assuming that 
the exact metric contains no high-frequency oscillations 
occurring on the body's natural timescale ~ e. In other 
words, the body is assumed to be in equilibrium. 

I seek a solution in a vacuum region outside the body. 
Given my assumptions, the vacuum EFE G = can be 
expanded as 

= G = Gilge] + 5Gi[H] + 6^Gi[H] + (56) 
where each term is further expanded as 

Gi[gB] = [d^\gB]+ eG^p[gB]+ e^Gf\gB]) , 

(57) 

6''Gi[gB] = e-^(^S''G^"\H]+eS''G^p[H]+eH''Gf'>[H]^ . 

(58) 

The overall factors of result from R = R/e and the 
fact that the Einstein tensor scales as the metric divided 
by two powers of length. The correction terms contain 
derivatives with respect to T, which arc each suppressed 
by a factor of e; specifically, G^"'' and S'^g\"'^ consist of 
the terms in G/ and S'^Gi that contain n derivatives with 
respect to T. Now, suppose H possesses an expansion 

Ni 

H{T,R,&^,e) = ^e"ij(")(T,^,e^). (59) 

n=l 

Substituting this expansion of H into the above expan- 
sion of the EFE, and then solving order-by-order in pow- 
ers of £, leads to the sequence 

Gf'"'[5i3]=0, (60) 
^^(o),,.j^(i)j ^ ^G^p^^'^[gg]^ (61) 

-Gf^^'iss], (62) 



Note that there is only one timescale here, so these equa- 
tions automatically follow from the assumed form of the 
expansion of the metric; there is none of the potential 
failings of a two-timescale expansion. In Sec. IV B and 
Appendix F, I discuss a particular solution to this se- 
quence of equations. 



IV. CALCULATION OF THE SELF-FORCE 
FROM MATCHED ASYMPTOTIC EXPANSIONS 

In this section, I consider the most intuitive means of 
solving the sequences of equations just presented: the 
method of matched asymptotic expansions. As outlined 
in the previous section, in this method the perturbation 
equations in the inner and outer expansions are solved 
independently, and then any free functions arc identi- 
fied by insisting that the two metrics agree in the buffer 
region around the body. Following the tradition of the 
field, in matching the two metrics I make use of the weak 
matching condition, rather than the strong condition. 

My presentation of the matching procedure roughly 
follows that of Refs. [37, 54], though most of my conclu- 
sions apply as well to the earlier calculation performed 
by Mino, Sasaki, and Tanaka [25]. However, my goal 
is not simply to review those earlier calculations, but 
to pinpoint their underlying assumptions. First among 
these assumptions is a very strong restriction on the re- 
lationship between the inner and outer solutions: essen- 
tially, the two solutions must be assumed to differ only 
by generically "small" coordinate transformations in the 
buffer region. This restriction is required because the 
weak matching condition, which has always been used in 
matched-expansion derivations of the self-force, is found 
to be too weak to yield unique results. The required 
restriction amounts to introducing a "refined" matching 
condition midway between the weak and strong condi- 
tions. 

Second among the underlying assumptions is the re- 
strictive choice of inner solution, which efltectively al- 
ready removes many of the integration constants that 
would normally be fixed by a matching procedure. As 
discussed in Sec. II, in traditional matched asymptotic 
expansions the leading-order inner and outer solutions 
are determined entirely by boundary conditions, while 
in the matched expansions used in the self-force prob- 
lem, the leading-order solutions must be chosen based on 
some desired physical properties; only after the leading- 
order solutions are chosen can boundary conditions be 
imposed. In the self-force problem, the leading-order 
outer solution is taken to be some desired vacuum met- 
ric. For the EMRI problem, the desired metric is that of 
a Kerr black hole. Typically, for simplicity, the leading- 
order inner solution is taken to be that of a Schwarzschild 
black hole, though one could instead choose, for example, 
that of a neutron star [57, 58]. 

However, in derivations of the self- force, the inner and 
outer solutions have been even further restricted: the 
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form of the perturbations have also been largely selected, 
rather than determined by matching. For example, the 
inner perturbations have been taken to be of a particu- 
lar form presumed to correspond to the influence of tidal 
fields on the small black hole. And the outer perturba- 
tion has been taken to be that of a point particle. In this 
section, I will make use of these assumed forms for the 
inner and outer solutions. As shown in Rcfs. [17, 27, 28], 
the assumption of a point particle perturbation can be 
removed, because the point particle solution follows di- 
rectly from the assumed existence of an inner expansion. 
On the other hand, the assumed form of the tidally per- 
turbed black hole metric has not, to my knowledge, been 
rigorously justified. Instead, I will point out the ways 
in which this metric restricts the generality of the inner 
solution. 

My analysis begins with a discussion of the outer ex- 
pansion. Section IV B then describes the metric in the in- 
ner expansion. In Sees. IVC-IVD, I perform the match- 
ing procedure, focusing on the restrictions that must be 
imposed to yield a unique result. I conclude the section 
with a discussion of the method. 

The calculations in this section make use of numerous 
computational techniques that are standard in the liter- 
ature: near-coincidence expansions, Fermi and retarded 
coordinate systems, and curved-spacetime Green's func- 
tions, STF (symmetric trace-free) decompositions, and 
tensor-harmonic expansions. Refer to Ref. [54] for ped- 
agogical reviews of the first three techniques. STF de- 
compositions and tensor-harmonic expansions are briefly 
outlined in Appendices D and F. 



A. Outer expansion 

I require an expansion of the background metric g and 
the first-order external perturbation /i^^ in the buff'er re- 
gion. To find these expansions, I adopt Fermi coordinates 
{t,x°') centered on 7 and then expand in powers of the 
geodesic distance r = y/Sijx'xi . The construction of the 
coordinate system is derived in Ref. [54]. 

Uppercase Latin indices /, J, K run from to 3, corre- 
sponding to an ortho normal tetrad e" . Lowercase Latin 
indices i, j, fc run from 1 to 3 and correspond to the spa- 
tial part of either the coordinate basis or the orthonor- 
mal basis. I am interested only in components in the 
Cartesian-type coordinates (t,a;°), but I will sometimes 
express these components in terms of r and two angles 
9^, which arc defined in the usual way in terms of x". 
I also introduce the unit one-form Ua = dar, which de- 
pends only on the angles 6^, and use the multi- index no- 
tation = n*i...n'* = n'^ -*'^. Angular brackets around 
indices denote the STF combination of the enclosed in- 
dices; a caret over a tensor denotes the STF part of 
that tensor. Finally, I define the coordinate one-forms 
ta = dat and X° = daX'^. 

In Fermi coordinates, the components of the back- 



ground metric are given by 

gtt = -1 - 2ra,n^ - ir^aiO* - a(iaj)n'^ 

-r^E^jh'^ +0{r^), (63) 
gta = yea^kB^n'^ +0{r^), (64) 

gab = Sab - ^r'^Sab£ijn'-^ " ^r'^^ab 

+ y£,^anl)+0{r^), (65) 

where I have decomposed the components into irreducible 
STF pieces, and defined the tidal fields £ab = RaObo and 
Bab = ^^a'^'^Robcd- Thc tidal fields are functions on the 
worldline, and are therefore functions of t only. 

One should note that the coordinate transformation 
x"(t, x") between Fermi coordinates and thc global coor- 
dinates is e-dependent, since Fermi coordinates are teth- 
ered to an e-dependent worldline. If one were using a 
regular expansion, then this coordinate transformation 
would devolve into a background coordinate transforma- 
tion to a Fermi coordinate system centered on a geodesic 
worldline, combined with a gauge transformation to ac- 
count for the e-dependence. But in the present general 
expansion, the transformation is purely a background 
transformation, because the e-dependence in the trans- 
formation is reducible to the e-dependence in the fixed 
worldline. 

The transformation hence induces not only new e- 

(n) 

dependence into the perturbations , but also e- 
dependence in the background metric g. (Despite its 
e-dependence, g is the background metric of the outer 
expansion, and I will use it to raise and lower indices 
on h.) This new e-dependence takes two forms: a func- 
tional dependence on 2"(<) = x°'{t,x'^ = 0), the coor- 
dinate form of the worldline written in the global coor- 
dinates x"; and a dependence on the acceleration vec- 
tor a"(t) on that worldline. For example, the first type 
of dependence appears in the components of the Rie- 
mann tensor (or tidal fields) in Fermi coordinates, which 
are related to the components in the global coordinates 
via the relationship RijKL{t) = Rai3-fs{z^{t))efeje'J^ej^. 
The second type of e-dependence consists of factors of 
the acceleration a'^{t), which has the assumed expansion 

a,(i,e)=af (i) + eaf)(t;7) + 0(e'). 

Hence, in the buffer region we can opt to work with the 
quantities g and He, which arc defined with a fixed, or 
wc can opt to re-expand these quantities by substituting 
into them the expansion of a. (In either case, we would 
still hold fixed the functional dependence on z^.) Sub- 
stituting the expansion of a in Fermi coordinates yields 
the buffer-region expansions 

Of.. = 5^.°^ {t, x''; 7) + eg^^Jit, x^ 7) + 0{e^) , (66) 
/^L1, = C^(^,^";7) + 0(e), (67) 

where gl°J = ff^^L^^m, gi^J = af^^f L=a(o)' and 
h\j^J = /i^'' I _ (0) • Because the inner expansion does 
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not hold the acceleration fixed, for the sake of matching, 
I will use the buffer-region quantities g'^^\ g^^\ and h^^\ 
In order to determine h^^\ I rewrite the wave equation 
(45) as 



-L II. 



(68) 



where Tajj is the stress-energy tensor of a point particle, 
given by 



■5'^{x~ z{t))dt, 



(69) 



wher 



dt 



is the four- velocity on 7, and \g\ de- 



notes the absolute value of the determinant of g^p ■ Note 
that there is no contradiction between this equation and 
Eq. (45), since the latter applies only in the vacuum re- 
gion 51, where Tap vanishes pointwise. The solution to 
this wave equation can be expressed in terms of an inte- 
gral over the worldline 7. Near the worldline, the solution 
can then be expanded in powers of r. That calculation 
is presented in Appendix E. The result is the following: 
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(70) 
(71) 



(1.0) 



.J, /(l.l), 
-i {a-^b)j ' 



(72) 



Here the uppercase hatted quantities are STF Cartesian 
tensors that are functions of time alone; they are named 
following the scheme of Eqs. (D14)-(D16). They are con- 
structed from tail integrals, the acceleration, and £, and 
their exact form is specified in Table I. Together, they 
make up the Detweiler- Whiting regular field [59], a so- 
lution to the homogeneous linearized Einstein equation. 
/j(i) is given by setting a,; = af''^ in the above expressions. 

The buffer-region expansion of the full metric in the 
outer limit can now be written as 

gE.p = gTp + + + 0(t\ er\ e^) , (73) 

where g^p is given in Fermi coordinates in Eqs. (63)-(65). 

B. Inner expansion 

I assume that the internal solution is that of a per- 
turbed Schwarzschild black hole, and I adopt retarded 



Eddington-Finkelstein coordinates ([/, X°) adapted to 
that spacetime. The background metric g^ is then given 

by 



gs ^ -f{U, R)dUdU - 2fladUdX'' 

+ {Sab - flab)dX''dX\ 



(74) 



where /(f7,^) = and = X'^/i? is a function 

of two angles 8^. (Note that in this equation, I have 
written the metric in non-rescaled coordinates, but I have 
written the components of the metric in terms of the 
scaled coordinate R.) Here M{U) is the Bondi mass of 
the spacetime, divided by the mass at U = 0. The mass 
is allowed to depend on U because gs is required only to 
solve Eq. (60), which contains no time-derivatives. Next, 
I expand the components of the metric perturbation H 
as 



H^^iU, R, 6^, e) = ei/W ^u, i?, 9^ 



) 



■'Hj^JiU,R,e^) + ... 



(75) 



As a boundary condition on these perturbations, I re- 
quire that they remain regular on the event horizon. In 
addition, I adopt the light cone gauge [60] , defined in re- 
tarded polar coordinates by the condition Hjj'^ = -ff = 



(») 

RA 



0. In this gauge, U and R maintain their geo- 



H 

metrical meaning even in the perturbed spacetime: U is 
constant on each outgoing light cone, and R is an affine 
parameter on outgoing light rays. I assume that this 
gauge condition can always be imposed. 

The first- and second-order perturbations, along with 
the time-dependence of gs, must satisfy the vacuum 
Einstein equations (61)~(62). In Appendix F, I show 
that ^ ~ 0. This implies that Eq. (61) becomes 

SC'-plH^^^] = 0, the linearized vacuum EFE for static 
perturbations. The solutions to this equation have been 
thoroughly studied [61-64]. Because of the spherical 
symmetry of the background spacetime, the equation can 
be most easily solved by decomposing iJ^^-* into spheri- 
cal harmonics: the various harmonics decouple in the 
linearized Ricci tensor, such that they can be solved in- 
dependently. In addition, for £ > 0, the harmonics can 
be decomposed into even- and odd-parity sectors, which 
also decouple. It is known that the gauge-invariant con- 
tent of the monopole terms in the solution corresponds to 
a constant shift of the black hole's mass parameter; odd- 
parity dipole terms correspond to a shift to a nonzero, 
constant spin; and even-parity dipole perturbations cor- 
respond to a shift in center of mass, which can always be 
removed via a coordinate transformation. In addition, it 
is known that for all £, the solutions behave as ^ R^ for 
^> 1. 

In the derivation provided by Mino, Sasaki, and 
Tanaka [25], the monopole and dipole terms in i/*-^' were 
set to zero, on the basis that they correspond to either 
pure gauge or to mere redefinitions of mass and angu- 
lar momentum. However, this step is not justified, since 
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the "constant" shifts in the black hole's parameters are 
actually functions of time, with a time-dependence to 
be determined by the higher-order perturbation equa- 
tions. Also, the fact that the even-parity dipole term 
corresponds to a shift in center of mass does not mean 
that it can be trivially ignored; this will be discussed fur- 
ther in the following sections. For £ > 1, Mino, Sasaki, 
and Tanaka took the terms to necessarily behave as in 
the buffer region. This means that H^"^ cannot contain 
terms of ^ > n: since e"i?^ = e"~^i?^, if ^ > n then such 
a term would correspond to negative powers of e in the 
outer expansion. Hence, H^^^^ can contain only monopole 
and dipole terms, and since these are set to zero, H^-^^ 
itself must be zero. It then follows that Eq. (62) becomes 
another linearized vacuum EFE for static perturbations, 
SCf''^ [H^'^'^] = 0. The solutions to this equation must 
be purely quadrupolar, since monopole and dipole terms 
are set to vanish and ff*^^^ cannot contain terms of ^ > 2. 
However, even if the monopole and dipole terms are set to 
zero, this reasoning remains specious, because solutions 
with £ > n can exist: though the asymptotically domi- 
nant terms behave as R^, subdominant terms can grow 
less rapidly with as is shown explicitly in Appendix F. 
(However, these subdominant terms can be removed with 
a gauge transformation.) 

In the derivation provided by Poisson Ref . [37, 54] , all 
of the above steps were taken, but the quadrupole terms 
were then further constrained. Rather than finding a gen- 
eral inner solution and then restricting it by imposing a 
matching condition, Poisson simplified the possible forms 
of the metric by first imposing a form of the asymptotic 
matching condition. (Refer back to Sec. II for the defini- 
tion of this condition.) Specifically, he demanded that for 
i? 3> 1, the metric must asymptotically approach that of 
a vacuum spacetime in retarded coordinates centered on 
a geodesic. This demand motivated the following ansatz 
for the internal metric in polar coordinates: 



giuu = -/ 
giuA = R 



l+e^R^ei{R)£*{U) 



O 



2^2 
3 



R^(^e2{R)£X + b2{R)B*x) + 0{ 



giRM = giRA = 0, 

giAB = i?' \^AB - k'i?' ( ei{R)£\B + h{R)B*AB 



(76) 
(77) 
(78) 
(79) 

) 

(80) 



of R could also be functions of C/, with a [/-dependence 
to be determined by the higher-order EFE. 

One might wonder why the internal metric is con- 
strained to approach that of a vacuum metric in coor- 
dinates centered on a geodesic, rather than being con- 
strained to approach a vacuum metric in coordinates cen- 
tered on an arbitrarily accelerating worldline, to agree 
with the form of the external background metric in re- 
tarded coordinates. The reason is that the terms linear 
in the acceleration in the external background metric are 
even-parity dipole terms, which have been set to zero 
to ensure that the coordinates are mass-centered. I will 
return to the relevance of these terms in the following sec- 
tions, but I will note now that this assumed form already 
suggests that the body must be moving on a geodesic 
of some spacetime. That spacetime will turn out to be 
g + h^, rather than g. See Ref. [65] for further discussion 
of this point. 

Substituting the ansatz into the linearized EFE and 
imposing regularity at the event horizon determines the 
free functions. After transforming the resulting metric 
back into Cartesian-type coordinates, one finds 



giuu = -f-£'.fR'£* + 0{e'), 
giua = -^a + h^R'f{£: + K) + 0{e'), 



1 ^2 



e'R' 1 



3" 

giab = Sab — ^ab 3 

- '^s^R'B:, + 0{e% 



''ah 



(81) 
(82) 

(83) 



B\^^, £. 



ab 



£ 



AB^a^b, and 



where £, 

^ab ~ ^AB^a^b- ^otc that there is no a priori relation- 
ship between the mass eM of the internal spacetime and 
the mass em of the point particle perturbation in the ex- 
ternal spacetime. Similarly, although the inner solution 
was specifically constructed to asymptotically approach 
the form of an external metric in the buffer region, there 
is no priori relationship between £ab and £ab or between 
Bab and Bab- These relationships are to be determined 
in the matching procedure. 



where ei, 62, 63, 62, and 63 are undetermined functions 
constrained to approach 1 for i? 3> 1, and the quantities 
£*, B\, etc., are constructed from tidal fields £ab and 
Bab, as displayed in Eqs. (F104)-(F107). At ^ ^ 00 (or 
£ = 0), this metric is precisely the metric of a vacuum 
spacetime in retarded coordinates centered on a geodesic. 
Note that even with the constraint on the asymptotic 
behavior of the internal metric, the above ansatz is more 
restrictive than it need be: generally, the free functions 



To expand the metric in the buffer region, we rewrite 
R as R/e and then re-expand in powers of e; this cor- 
responds to an expansion for R ^ e. In order to agree 
with the external metric, which is constructed in Fermi 
coordinates and in the Lorenz gauge, we must also trans- 
form from retarded coordinates and the lightcone gauge 
into Fermi-like harmonic coordinates (T, X°')] and the re- 
sult must be decomposed into its irreducible STF pieces. 
That calculation is shown in Appendix F. The final result 
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(84) 



giTa = 2eMREa,N' + ^eM Rea^JBiN'^ 



(85) 
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2M 



ab ^Sab{l + e— - hMR£,,N'^ - \R^£,.N'^ 



+ ^eMRE^^Nt)' - §eMR£ab - ^R^Eab 

+ leMRS.jNab'' + lR%{a% 

- \eMRe^k(a%N^ + 0{e^,eR^,R') , 



(86) 



where iV* — X'^jR. This is the metric that I will use 
in the matching procedure, even though, as pointed out 
above, it has already been heavily restricted. 

In Appendix F, I derive the above metric without re- 
lying on an ansatz, enabling me to better characterize its 
generality. 



C. Zeroth-order matching 

I now consider the relationship between the two met- 
rics. Beginning with the zeroth-order weak matching 
condition, we have the metric in the outer expansion 
given by 



(0) i 12 (0) (0)1 2 (0) (0)~i,- 
3 J (i j> 



(87) 



gSab = Sab - Ir^SabStjfl'^ " ^r^^ab + ^r'^Si^a'^l^ 

+ 0{e,r^). (89) 
while the metric in the inner expansion is given by 



giTT = -l- R^E^jN'^ + 0(£, R"), 



giTa = ^R^earkB'^N'' + 0{e,R'') 



(90) 

3-^^- ^aiti'^j^-' ' y, (91) 

giab = Sab — \R^Sab£ijN^^ — \R^£ab + ^R^^iia^b) 

+ 0{e,R^). (92) 

It seems that we may immediately identify these two met- 
rics and conclude that T = t, X"' = x"', Sab = Sab, 
and most importantly, a|i*''' ~ 0. However, the match- 
ing condition does not require that these two metrics be 
identical, since they may be in different coordinate sys- 
tems; the matching condition requires only that these 
two metrics be related by a diffeomorphism. But this 
condition places no restriction at all on the acceleration 
of the worldline: The form of the inner metric is that 
of an arbitrary background written in Fermi coordinates 
centered on a geodesic worldline. The form of the outer 



metric is that of a known background written in Fermi 
coordinates centered on a possibly accelerated worldline. 
Regardless of the value of the acceleration, if the geodesic 
is embedded in the external spacetime, then these two so- 
lutions are obviously related by a diffeomorphism, since 
the geodesic can be transformed to the accelerated world- 
line. 

Evidently, some information has been lost here. I as- 
sumed from the beginning that the inner and outer ex- 
pansions were performed "around" the same worldline. 
In the inner expansion, the "location" of the body is en- 
coded into the coordinate system by the condition that 
the body's mass dipole vanishes in that coordinate sys- 
tem; in the outer expansion, the "location" of the body 
is encoded in the worldline sourcing the perturbation. If 
we use the weak matching condition, in which we expand 
the metric before finding the coordinate transformation 
between the inner and outer expansions, then this infor- 
mation is lost. 

However, one might wonder if this ambiguity might be 
removed by supplementing the weak matching condition 
with some other condition. One such condition appears 
obvious: the coordinate transformation between the in- 
ner and outer expansions in the buffer region must be 
"small" -that is, it must vanish in the limit e — > 0. This 
removes the possibility of transforming from an arbitrary 
geodesic to an arbitrarily-accelerated worldline. In the 
buffer region, r — as e — > 0, so this allows transfor- 
mations that have no explicit e-dependence, but which 
do have explicit r-dependencc. I trust the reader to con- 
vince himself that under such a transformation, we must 
have R = r, T = t, the tidal fields appearing in the in- 
ner metric must be identical (up to 0{s) corrections) to 
those constructed from the Riemann tensor in the outer 
solution — and the leading-order term in the acceleration 
must vanish: a'-^' = 0. The two coordinate systems may, 
of course, be related by rotations, but these are insignif- 
icant. 

Hence, we can adopt a stronger, refined matching con- 
dition: the inner and outer expansions in the buffer re- 
gion must be equal up to a unique small coordinate trans- 
formation. Unfortunately, this refined condition is still 
insufficient. The reason is that the inner expansion could 
have included acceleration-type terms. In fact, we can 
always include such terms by transforming the metric 
into an accelerating frame. Suppose we begin with the 
Schwarzschild metric in Kerr-Schild form. 



eM 



(93) 



where 77 = diag(— 1, 1, 1, 1) is the Minkowski metric, 

is a nuU vector, i? = \/X^ + + Z^, 

and the (unsealed) coordinates are {T,X,Y,Z). Now, 
by using the flat-spacetime transformation from an iner- 
tial frame to an accelerated one, we can transform the 
metric to a new set of accelerated retarded coordinates 
{U',R',Q'^). For simplicity, assume that the accelera- 
tion is in the Z-dircction. Then the transformation is 
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given by 

T = Tq{U') + R'icose' smhq{U') + coshq{U')), (94) 
X = i?' sine' cos (95) 
r i?' sine' sin (96) 

Z = Zo(C/') + i?'(cose'cosh(7([/') +sinhg([/')), (97) 

where To = / coshq{U')dU', Zq = J smhq{U')dU', and 
q{U') ~ J a{U')dU', where a{U') is the magnitude of the 
acceleration. Under this transformation, gB maintains 
the form in Eq. (93). ry^^ becomes the metric of flat 
spacetime in retarded coordinates, given by 

r]^-[{l + R'a cos Q'f - R'^a^] dU'^ - 2dU'dR' 
- 2R'^a sin O'dU'dS' + R'^dfl'^, (98) 

while takes on a more complicated (and unenlight- 
ening) form. Note that in flat spacetime, this trans- 
formation translates the spatial origin from Z = to 
Z = Zq{U'). And in the spacetime oi gs, the same inter- 
pretation applies at large distances from the black hole — 
that is, in the buffer region. In other words, the new 
coordinates arc not mass-centered: the center of mass is 
moving away from the origin. 

Although the metric takes on an inconveniently com- 
plicated form in this non-mass-centered coordinate sys- 
tem, in principle one could use it in constructing the inner 
expansion gj. If one does so, then when gi is expanded 
in the buffer region, it becomes Qi^y — 
as we can infer immediately from the form of Eq. (93). 
But in this expansion, rj^^ is the metric of flat spacetime 
centered on an accelerating worldline, not on a geodesic. 
Therefore, if we transform the metric to Fermi-type co- 



ordinates (T, X"), we arrive at 

5/TT = -1 -2i?a cose + 0(£,i?2), (99) 

5/Ta = 0(e,i?2), (^00) 

giab^5ab + 0{e,R^). (101) 



This metric agrees with the one in the outer expansion, 
regardless of the value of the acceleration. We may iden- 
tify a cose with af^n^, and the matching procedure, 
even with the refined matching condition, provides no 
information whatsoever about the worldline. 

We can readily see why the matching procedure has 
failed: we have not insisted on any relationship between 
the inner and outer expansions. In order for matching to 
be successful, we must insist that the "position" of the 
black hole in the inner expansion can be identified with 
the position of the worldline in the outer expansion. To 
make this identification mathematically precise, I insist 
that the two expansions are to be expanded and matched 
in the buffer region only when the outer expansion is eval- 
uated in a coordinate system centered on the worldline 
and the inner expansion is evaluated in a mass-centered 
coordinate system. If this condition is imposed, then the 
accelerating coordinate system (J7', i?', e^^) is inadmis- 
sible, since it is not mass-centered. Therefore, we can 



discount it and others like it — and we can once again, 
now more confidently, conclude that the acceleration of 
the worldline must vanish in the limit £ — > 0. Such a con- 
dition serves to implicitly define the worldline, and it is 
necessary for the matching procedure to be well-defined 
and to yield unambiguous results. 

Based on the above analysis of the zeroth-order match- 
ing procedure, I suggest the following matching condi- 
tion: if the inner expansion is written in a mass-centered 
coordinate system and the outer expansion is written in 
a worldline-centered coordinate system, then the two ex- 
pansions must be equal up to a small coordinate trans- 
formation when expanded in the limit of small (outer) 
radial coordinate distances. (Here "outer" radial coordi- 
nate means a coordinate that is formally of order 1 in the 
outer expansion and of order l/£ in the inner expansion.) 
Making use of this condition allows us to determine the 
acceleration of the worldline at zeroth order. However, as 
we shall see in the next subsection, it requires still more 
restrictions. 



D. First-order matching 

Comparing the expression for the external solution 
with that for the internal solution, we find that the 1/r 
terms agree if and only if we make the identification 
m = M. In order for the other terms to be made to 
agree, there must exist a coordinate transformation, from 
the external coordinates to the internal coordinates, that 
induces a gauge transformation gE gE + sSgs + Oi^e'^), 
where 

SgEu - -A^'^"^ + r(2aW - + 0{r^) , (102) 

SgEta = ~Ci'-'^ - ira,(i(i-°) + 37^(i'«))n, - rCj^'^^n^ 
~ rtoA't^^v} + 2mr£a^n' + 0{r') , (103) 

-I- ^mrSab + ImrSijfiab' 

- |mre,fc(,e^)n^ + 0{r^) . (104) 

I remind the reader that is to be set to af^ = in the 
explicit expressions for the uppercase script tensors. 

This transformation is generated by a vector field sat- 
isfying SgEap — 2^(„;^). I assume the field can be ex- 
panded as 

6 = E„>o^"d"'. ia^En>or''d"\ (105) 
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where the coefhcients ^j"-* and (,a"' arc decomposed as 
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the second of them yields 
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5^(1.1) 
8^ab 



The Cartesian tensors S^, T^, A^. and are STF in 
L, and they depend only on time. 

Calculating 2^(q.^) from the above expansion is 
straightforward. Demanding that the result of this cal- 
culation agrees with Eqs. (102)-(104) at each order in 
r then determines a sequence of equations for No 
0(l/r) terms appear in Eqs. (102)-(104), so from the 

.(0) 



+ |6^(,6^^Tf, (116) 
a,A(l) = -^(I'l) + i6/'?6',,6;^Tf\ (117) 
and the last of them yields (after some algebra) 



0{\/r) terms in 2£,^a.p) we find that daCt'' = and T^^^,^ = -Ti??'^^ 
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-Y-(2) 
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^ j ^ 16^a 8'^*^a ' 



= 0. From this we determine that £}a must be 

independent of angle: sjf^ = S^^^ and = Ti"-*. 
From the 0(r*') terms, we find 
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atsw = (108) 

(na+r9a)d''=-5*d"^-C(i'"\ (109) 
("(a + ra(j4^' = -i<5.bA'(i^") - (110) 

from the tt-, ia-, and a6-component, respectively. 
The first of these equations determines that S'"^ = 
— i /A'-^'^-'di, the second determines that si^-* = 



(1,0) 



9tTi°\ and the last determines that 4''^^' 



T^^^ = and Ai^^ is arbitrary All 

other terms in £}a ^ vanish. 

Finally, from the 0{r) terms, we find: 



(1) 



(111) 



(0) 



(112) 



2(2n(, + ra(JClf = ti?o(a6),n'S(o) - |i?^(,6),n^Tf ^ 



n' + ^mEii^a'hb)'' 



7 

ahi 



(113) 



^(2) 



4 ahc 



-iL(A:^i) + iiii^i)-ia,cr))-fif^T 



(2) 
(2) 



(118) 
(119) 
(120) 

(0) 

i ' 
(121) 

(122) 



if(i,i) 

2 ab 3 



2^ J 



(123) 



All other terms vanish. 

In summary, the first three terms in the expansion of 
the gauge vector field are given by 
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where T^"^ is a function of time satisfying the equation 
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of motion (114) 
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and 

f(2) 



ia,i^i^°) + 2£„S(0)+C^,^^^+2m£„ 
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1^(1,1) 



lij(i^i), 



+ \if,''> + ImB,, \B,u [A^'''^dt) . (129) 



Again, these equations follow from the tt-, ta-, and ab- 
component, respectively. The first of them yields the 



This is the most general transformation that succeeds 
in making the exterior solution identical to the interior 
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solution, up to order er. It has one free function of time: 

^ a ■ 

Despite the refinement of the matching condition for- 
mulated in the zeroth-order matching procedure, this co- 
ordinate transformation has failed to uniquely identify 
the acceleration of the worldline. Instead, it determines 
an equation for Ti°\ given by Eq. (114). Consider the 
meaning of this equation. In the internal solution, the 
mass dipole and all dipole perturbations have been set 
to zero, and an acceleration term in the buffer region 
corresponds to a dipole perturbation. Equation (114) 

thus tells us that for any given acceleration we can 
perform a small, angle- and r-independent spatial trans- 
lation (in the buffer region) that ensures the dipole per- 
turbation vanishes. 

In order to arrive at the correct equation, a'p = 
i^l^'") — dtC[^''^\ one must further restrict the match- 
ing condition. Recall that the coordinate transformation 
must be small, which implies that T^f^ must remain of 
order unity. If the right-hand side of Eq. (114) does not 

vanish, then t'"' will generically grow large; more pre- 
cisely, on a timescale such as ~ 1/e, which becomes un- 
bounded in the limit e — >■ 0, t|°-' will generically become 
larger than order unity. However, it will not necessar- 
ily grow large (e.g., if the right-hand side of Eq. (114) is 
purely oscillatory). Thus, we cannot conclude that the 
right-hand side must vanish based on the refined match- 
ing condition of the previous section. Instead, I propose 
a final version of the refined matching condition: if the 
inner expansion is written in a mass-centered coordinate 
system and the outer expansion is written in a worldline- 
centered coordinate system, then the two expansions must 
he equal up to a necessarily small coordinate transforma- 
tion when the inner expansion is expanded in the limit of 
small ( outer) radial coordinate distances. In other words, 
the coordinate transformation must not only be small, 
but must necessarily remain so on all timescales of inter- 
est. 

With this final refinement, we can conclude the follow- 
ing: if on an unbounded timescale, (i) the exact metric 
possesses inner and outer expansions, (ii) there exists a 
local coordinate system in which the metric in the inner 
expansion is given by that of a tidally perturbed black 
hole, up to errors of order e^, (iii) there exists a global 
coordinate system in which the metric in the outer ex- 
pansion is that of the external background g plus a point- 
particle solution to the wave equation (45), up to errors of 
order e^, and (iii) the exact solution satisfies the refined 
matching condition presented above, then the worldline 
defining the point particle perturbation has an accelera- 
tion given by 

„W^i^(M)_ 5^^(1.0)^ (130) 

where and cf"'^^ are obtained by setting Oi = 

af ^ = in Table I. This is the MiSaTaQuWa equation. 
It may be that some or all of these assumptions can be 



removed, even within the context of matched asymptotic 
expansions. For example, if the inner and outer expan- 
sions exist, then the solution to the wave equation (45) 
is necessarily that with a point particle source. Also, the 
inner metric could correspond to a body other than a 
black hole. If it were taken to be a tidally perturbed, 
otherwise spherically symmetric neutron star, for exam- 
ple, then the equation of motion would be unaffected: 
outside the star, the metric would be altered only by the 
presence of induced tidal moments, which scale as e^^"*"^ 
and hence would not appear in the first-order matching 
procedure [57, 58]. It is also possible that matching the 
inner and outer solutions to higher order in r would show 
that T^a^ must vanish, and that the refined matching con- 
dition is needlessly strong; however, there is no obvious 
indication of the order at which this would occur. 



E. Interpretation and commentary 

Let us interpret the above calculation. First, note 
that a large part of the transformation consists of re- 
moving tail terms. This can be understood as follows: 
In the Fermi coordinates centered on the worldline in 
the spacetime appears to be that of a singular monopole 
perturbation h^ , plus a regular homogeneous perturba- 
tion h^\ plus the field of the smooth background metric 
g expanded about some worldline. But in the coordi- 
nates X", at a large distance R ^ e from the body, 
the spacetime appears to be simply a singular monopole 
perturbation atop some smooth background field that is 
expanded around a geodesic. Therefore, transforming 
between these coordinates can be understood as trans- 
forming from the Fermi coordinates of g into the Fermi 
coordinates of g-\-h^, where is the Detweiler- Whiting 
regular field. Reference [65] contains further discussion 
of this point. 

For example, the angle- and r-independent monopole 
term A^^'^^ = hlf^^ is removed because proper time must 
be measured in g + , rather than in g. Similarly, the 
dipole terms in the perturbation are removed because 
the body is nonspinning and non-accelerating mg + h^, 
rather than in g. And if we proceeded to order er^ in the 
matching procedure, we would find that the tidal fields 
appearing in the inner expansion are those of g -\- h^, 
rather than those of g. 

These general points are shared with the derivation in 
Refs. [37, 54]. Note, however, that my results differ con- 
siderably from those of Refs. [25, 37, 54]. The first differ- 
ence is that in those earlier calculations it was found that 
the tetrad on the worldline is not parallel-propagated in 
the external background spacetime. This followed from 
the fact that the spin dipole perturbation in the internal 
solution had been set to zero via a choice of gauge; effec- 
tively, the Fermi tetrad was required to rotate with the 
perturbed gravitational field, to set the total spin to zero. 
However, this is not necessary: transforming from the ex- 
ternal, nonrotating frame, to the internal, rotating frame 
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(or to a nonrotating frame in g + /i^), simply requires 
a gauge transformation (specifically generated by Ac^'). 
There is no reason to require the external, background 
Fermi tetrad in g to spin. 

More importantly, my analysis has shown that the 
weak matching condition that is normally utilized is actu- 
ally too weak to yield unique results. In order to arrive at 
an equation of motion, I have had to formulate a refined 
matching condition in a somewhat ad hoc manner. But 
even if that refinement is accepted, all of these derivations 
rely on another strong assumption: the form of the in- 
ner expansion, which fixes not only the background met- 
ric, but also the form the perturbations. By assuming 
that the metric in the inner expansion appears to that 
of a singular monopolc perturbation atop some smooth 
background field that is expanded around a geodesic, the 
matching derivations of the self-force seem to implicitly 
assume a generalized equivalence principle: they assume 
that in vacuum, the black hole, as viewed from a distance, 
moves on a geodesic of some smooth spacetime. Given 
that such a geodesic exists, the matching procedure pro- 
vides a means of determining which smooth spacetime 
the geodesic lies in — but it does not prove the existence 
of the geodesic. 

If the inner expansion is to be sufficiently general for 
the matching procedure to derive the generalized equiva- 
lence principle, rather than assume it, then one must use 
a less restricted inner expansion. For example, at linear 
order, no acceleration-like dipole term (i.e. one behaving 
as ~ r in the buffer region) can arise in the inner expan- 
sion without also introducing a mass dipole. However, 
in the inner expansion an acceleration term era^^^ cor- 
responds to a second-order perturbation e^i?a'-^-', so in 
order to maintain that no such term can arise in mass- 
centered coordinates, one must solve the second-order 
EFE in generality. In addition, if one uses the refined 
matching condition or some variant of it, one must prove 
that an acceleration term cannot arise from a generically 
small coordinate transformation in the buffer region. It 
is not immediately clear that at second order, being in 
mass-centered coordinates implies that any acceleration- 
like term must vanish.^ Solving the EFE, in full gener- 
ality, at second order also requires one to include po- 
tentially time-evolving shifts in the mass and spin of 
the black hole. The time-evolution, if any, of these pa- 
rameters would be determined at second order. Indeed, 
a time-dependent correction to the mass was found in 
Ref. [28]. 

In that paper, I performed a second-order analysis, but 
of a different sort than the one just suggested. Instead 
of assuming that the small body is a black hole, and try- 
ing to solve the second-order EFE in the whole of that 



Kinnersley's photon rocket is an example of an exact solution to 
the EFE with acceleration-like dipole terms sourced by radiation 
[66-68]. 



black hole's spacetime, I allowed the body to be arbi- 
trarily structured, and I solved the second-order EFE in 
the outer expansion, and only in the buffer region. (See 
Refs. [15, 27] for similar approaches at first order and in 
the case of a regular expansion, respectively.) Given the 
weaknesses of the method of matched asymptotic expan- 
sions, such a calculation provides a much firmer conclu- 
sion. However, if one is willing to accept the additional 
assumptions, matched asymptotic expansions provide a 
much simpler means of finding an equation of motion. 



V. CONCLUSION 

In this paper, I have discussed three types of general 
expansions: dual expansions, multiscale expansions, and 
a self-consistent expansion. Each of these can be used to 
systematically overcome the unbounded errors of a regu- 
lar expansion, and at least two of them will likely be re- 
quired in any successful analytical approach to the prob- 
lem of motion for asymptotically small bodies. For exam- 
ple, whether a multiscale expansion or a self-consistent 
expansion is used to eliminate secular errors, it must be 
combined with an inner expansion to eliminate errors 
near the small body. And these expansions may need 
to be supplemented with a third expansion in the wave 
zone. 

However, these expansions are not merely powerful 
tools for finding asymptotic solutions to differential equa- 
tions: they have a rich underlying geometrical structure, 
which I have sketched, but which warrants further study. 
This underlying structure allows us to, for example, eas- 
ily define and conceptualize the representative worldline 
of a black hole. Careful consideration of the underlying 
formalism can also reveal subtle yet important features 
of a calculation. For example, in moving from the realm 
of traditional perturbation theory to one with potentially 
multiple coordinate systems, two types of matching con- 
ditions arise, and the condition that has been implicitly 
used in previous derivations of the gravitational self- force 
is significantly weaker than the condition used in applied 
mathematics. In order to arrive at unique results with 
this matching condition, additional assumptions must be 
made, which weakens the conclusions of a matching cal- 
culation. 

It is my hope that this paper, by highlighting points 
such as these, will spur further research into the founda- 
tions of singular perturbation techniques in GR, in order 
to promulgate their utility, illuminate their underlying 
structure, and put them on a more precise and deductive 
mathematical basis. 



Acknowledgments 

Thanks to Eric Poisson and Bob Wald for helpful feed- 
back, and to the anonymous referee for a careful reading 
of the mathematical expressions. This work was sup- 



25 



ported by the Natural Sciences and Engineering Research 
Council of Canada. 



Appendix A: Illustrative example of matched 
asymptotic expansions 

To illustrate the procedure of matched asymptotic ex- 
pansions, I consider the following boundary value prob- 
lem: 



solving order-by-order, we find 



(A4) 



+ ^+f = 0, f(0) = 0, f(l) = l. (Al) 
ar 



We can now make use of one of the matching conditions 
to determine the integration constants Z)*-°^ and _D'^^ 

Since limf _>oo //„ ^ (r) = -D^") and lim,._^o /iut = e, the 
asymptotic matching condition implies Z?*-*^^ = e. In or- 
der to determine D'^^\ I next make use of the coefficient- 
matching condition. Rewriting as a function of r and 
expanding to order e, we find 



I assume that f possesses an outer expansion /out = 

/out('') + £/out(^) + Substituting this expansion into 
Eq. (Al) and equating coefficients of each power of e to 



/in(r) 



er 



(A5) 



zero, we find + f^^^ 



AO) 



and 



the boundary conditions are /out(O) 



dr 



J out 



/out(l) = 1; ^^'^ /out(l) = 0- The general solution to the 
zeroth-order differential equation is f^nti''^) = C^^^e"*". 
This solution can satisfy the boundary condition at r = 1 
(by setting C^^^ = e), but it cannot satisfy the condition 
at r = 0. Hence, we guess that there is a boundary layer 
at ?■ = 0, and we choose only to satisfy the boundary 
condition at r = 1. Doing the same for /quL find 



f(i) 



/i;i(o) 



0, 



(1) 

J out 

sion 



(1 — r)e^ ^, yielding the first-order outer expan- 



/out -ei-'^+e(l-r)ei 



(A2) 



Note that e^''^*^ — o(e") for all n > 0, so it vanishes 
in this expansion. Dropping the ellipses, this yields the 
left-hand side of Eq. (3): 

^ItP*<PItPJ ^ e + sD'^^^ ~er. (A6) 

Next, expanding /out to linear order in r, we find 

/out = e(l-r)+e(l-2r)e+... (A7) 

This expansion contains an order-£r term, which is 
smaller than any term in Eq. (A5); such a term would 

(2) 

be matched by a term from /;„ , and so we can neglect 
it here. The extra expansion on the right-hand side of 
Eq. (3) serves to remove such terms, and so we have 



$^i/'*$^i/'*$^f = e(l - r) + ee. 



(A8) 



Now, in order to construct our inner expansion, we re- 
quire a choice of rescaled coordinate r. Suppose that we 
choose r = r/e^. Substituting this into Eq. (Al) and 
taking the limit e — )■ 0, we find that if < p < 1, then 



Hence, the matching condition now determines that 
= e, and we have fully determined the inner ex- 
pansion: 



df'-° 

the leading-order differential equation is ^ = 0. If f^^ = e(l - e^'') + se [(1 - e"*^) - f(l + e"*^)] + ... (A9) 



p > 1, then the equation becomes = 0. And if 

p ~ 1, then it becomes — — I- = 0. This is called 
a distinguished limit (also known as a significant degen- 
eration) of the equation, because it contains within it all 
the terms appearing in the other two limiting equations. 
We can intuit that a distinguished limit will yield an ap- 
proximation with maximal information. And although 
there is no guarantee that a coordinate leading to a dis- 
tinguished limit is the ideal choice, it has proven to be 
the most reliable one. 

So, proceeding with the rescaled variable f = r/e, I 
rewrite Eq. (Al) as 



Using these results, we can construct the uniformly 
accurate composite expansion 

./comp = {e'-'- + £(1 - r)e^-^} + {e(l - e-^/') 
+ ee (1 - e-''/"') - r/e(l + e-"^/") } 
— {e(l — r) + ee} 



- (l + r)ei-'^/^ 



(1 



(AlO) 



rf2f 



df 



z^ + -Tz+ef = 0, f(0)=0. 



(A3) 



where the first equality should be compared to Eq. (6). 

In this case, we can compare our results to the exact 
solution to Eq. (Al), which is given by 



where, technically, f stands in for -0*f. I assume that 
f(f) possesses an inner expansion /in(r, e) = fi^\f) + 
ef^^\r) -I- .... After substituting this into Eq. (A3) and 



f = 



exp [-(1 - VT^)^] - exp [-(1 + VT^)^] 

exp [-(1 - VT^) Jj] - exp [-(1 + VT^)^] ' 

(All) 
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FIG. 5: Comparisons of the exact solution f (the sohd black 
curve), the inner solution /in (dotted blue), the outer solu- 
tion /out (dot-dashed red), and the composite solution /comp 
(dashed black). The upper plot displays the solutions for 
e = 0.2; the lower, for e — 0.1. 



Note that this function does not exist at e = 0. Hence, 
the regular series /out is not a Taylor series expansion of f . 
However, the limit limg^o f does exist, and /out is given 
by ^ |j limg^o fjl-- One can straightforwardly check 
that /in (when written as a function of r) is a uniform, 
first-order asymptotic approximation on an extended do- 
main Din = {r : < r <^ 1}, and /out is a uniform, first- 
order approximation on Dout = {r : |elne|<g;r<l}. 
So at first order, the overlap region exists, and it is given 
by I £ hie I <C r <C 1 (which is notably smaller than the 
buffer region). One can also verify that /comp is a uni- 
form approximation on the whole interval [0, 1]. Figure 5 
shows a graphical comparison of the exact, inner, outer, 
and composite solutions. 



Appendix B: Illustrative example of multiscale 
expansions 

In this appendix, I present an illustrative example of 
multiscale expansions, along with a resultant discussion 
of their utility. I consider the following differential equa- 
tion, adapted from the text by Kevorkian and Cole [44] 



f(0,£)=0, f (0,£)-l. (Bl) 



Suppose we wish to solve this problem using a regular 
power series f{t,e) = E„>o^"/^"H*)- After substituting 
this series and equating powers of e, we arrive at the 
sequence of equations 



d2/W 



+ /(0) 
+ /(!) 



0, 



dt '' 



/(o)(o) = o, ^(0) = 1, 
/W(0) = 0, ^(0) = 0. 



(B2) 

The solutions to these equations are easily found to be 
/(o)(t) — sint and f^^''{t) = — tsini, so we have 



f(i,£) 



sin t — et sin t 



(B3) 



Based on the unbounded growth of this solution, we sur- 
mise that it fails to uniformly approximate the exact so- 
lution on any unbounded interval [0, p > 0. 

To improve on this solution, I adopt the following as- 
sumption: there exists a function F(t, t, e) satisfying the 
equality F{t,i = et,e) = f{t,e). Substituting this into 
Eq. (Bl) and making use of the chain rule ^ = ■§t+£-§^, 
we arrive at 



fd'^F 


dF\ 


..( 


\dtdt ^ 







d^F 



dF 



dt 



+ F + 2e\ -:— + — I + e I — + 2— ) = 0. 

.(B4) 

Now, the fundamental idea in a multiscale expansion is 
that the function F satisfies this equation not just when 
t = et, but also when t is treated as an independent 
coordinate. This means that if one assumes a regular 
expansion F{t,i,e) = X]n>o ^"^^"''('•' then the coef- 
ficient of each power of s in Eq. (B4) must vanish; if 
we insisted on solving the equation only at t = et, then 
the e-dependence embedded in t would prevent us from 
concluding that the equation must be satisfied order-by- 
order in this way. (Of course, we could always solve the 
equation by setting the coefficient of each power of e to 
zero, but we could not deduce that each coefficient must 
vanish.) 

So, following this procedure, we arrive at a new se- 
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qucncG of equations, 

^2^(0) 



9t2 
52 ^(2) 



+ = 0, 

9^(0) Q2p(0) 

= -2— 2- 



9i 



ai^(i) 92j^(i) 92^(0) 
F^^^ = -2 2- 



9t 



9t2 



9t 



(B5) 
(B6) 

(B7) 



subject to the initial conditions i^^"^(0,0) = for n > 0, 

^(0,0) ^ 1, and ^(0,0) = -^^^(0,0) for n > 
0. The solution to the first equation is 

^(0) ^ ^(0) gij^ t + (t) cos i, (B8) 

where the initial conditions on F*^*^' do not fully deter- 
mine the slow evolution of A^^^ and B'-^\ but only impose 
A(°)(0) = 1 and B(o)(0) = 0. The general solution to the 
second equation is 

pW = A^^\i)smt + B^^'>{i)cost 

A(0)(i) + ^^-^(i) (cost + 2tsint) 
dt ) 

(f) + (t) t cos t. (B9) 

/ 

I now make a final assumption, called the no-secularity 
condition: the ratio of successive terms, F^'^~^^'> / F^^\ 
must be bounded. This means that the terms isint 
and t cos t are inadmissible, and so their coefficients must 



vanish. In this case, we have 



dt 



and 



dt 



+ ^ 0, subject to A(o)(0) = 1 and B^°\0) = 0. 

-* and = 0. 



Solving these equations, we find A^^^^ = 
Hence, we have now fully determine F^^' to be 



pio) 



e *sint. 



And we have 



^(1) = A^^\i) sin t + B^^\i) cost, 



(BIO) 



(Bll) 



where the initial conditions on F'^^^ imply that ^'^^(0) = 
B(i)(0) = 0. 

If we ceased our work here, there would be no signal 
that our assumed expansion cannot, in fact, satisfy the 
non-secularity condition. Following the same procedure 
for Eq. (B7) as we did for Eq. (B6), we find that in order 
to avoid secular growth in F'-^\ the functions A^^^ and 



B^^) must satisfy the equations 



^(1 



= 



and 



dB^ 



dt 



0; along with the initial condi- 



tions yl(i)(0) = B(i)(0) = 0. The solutions to these equa- 
tions are the secularly growing functions A^^-* = ie~* and 



5(1) 



Thus, in order to avoid secular growth 



m we must introduce secular growth into F^^K In 

other words, the expansion has failed. 

In this case, we can determine the precise reason for 
the failure. The exact solution to the original ODE is 



f(t,£) 



-et 



: sm 



(B12) 



If we expand this in a regular power series, we arrive at 
f(t,s) = sin< — eisini -I- agreeing with the regular 
expansion given in Eq. (B3). But we find by inspection 
that f (<, e) cannot be written as F(t, t, e) in such a way 
that a regular expansion of F satisfies the no-secularity 
condition. While e"*^* can be written as e^* to remove 
secular growth, an expansion of sin(t-\/l — s^) will violate 
the condition. However, we can write f(t,£) = F{(j),t,e), 
where (f) ~ ^{s)t, 51(e) = \/l — e^. and F is given by 



^(0,f,e) 



■ sm( 



(B13) 



This function possesses the regular expansion F{(j>, i, e) 



(1- 



' sin (f>+o{e), which, when expressed in terms of 



i, is a uniform approximation to f{t,e). One might won- 
der if we could have discovered this expansion without 
access to the exact solution. The answer, fortunately, is 
that we could have: substituting f = F = J2 £"F^"''> {(j), i) 
and n{e) = J2n>Q e"^^*"'' into Eq. (Bl) and then solving 
for arbitrary (j) and t yields a sequence of equations that 
determine the F^") and Q^"^ [44]. 

There are several points to note from this example. 
First, while an expansion method might appear to be 
working, it might still fail at higher order. Second, al- 
though we cannot be guaranteed that this failure will re- 
veal itself in the course of our perturbation calculation, 
that will typically be the case, as it was here. Third, 
even though my assumptions about F proved to be false, 
and even though F^*^) -|- eF^^^ fails to provide a uniform 
first-order approximation to f, the term F^^'f alone, the 
only term in F that was fully determined without any ob- 
vious contradiction, does provide a uniform zeroth-order 
approximation. (This can easily be checked by calculat- 
ing the supremum norm of |f — F*-*^'!.) 



Appendix C: Gauge transformations in the 
self-consistent expansion 

In the self-consistent expansion of Sec. HIE, the outer 
expansion is defined not only by holding fixed, but 
also by demanding that the mass dipole of the body van- 
ishes when calculated in coordinates centered on 7. If we 
perform a gauge transformation generated by a vector 
^W°'(^x;^), then the mass dipole will no longer vanish in 
those coordinates. Hence, a new worldline 7' must be 
constructed, such that in coordinates centered on that 
new worldline, the mass dipole vanishes. In other words, 
in the outer expansion we have the usual gauge freedom 
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of regular perturbation theory, so long as the worldline is 
appropriately transformed as well: (ft,, 7) -> (ft.', 7'). The 
transformation law for the worldline is well known [69] ; it 
will be worked out again presently. The only new feature 
of this gauge freedom is that e-dependence can be incor- 
porated into the transformation, because the gauge vec- 
tors can be functionals of the old worldline; this allows, 
for example, a "first-order" gauge vector that is con- 
structed from the tail integral ft*'^''[7]. However, to main- 
tain the form of the expansion, we must also insist that 
^(n) _ 0^(1). Of course, in addition to this gauge free- 
dom, one can still perform global, e-independent back- 
ground coordinate transformations. 

I first justify, to some extent, the assumption that the 
Lorenz gauge condition can be imposed on the entirety 
of ft. If we begin with the metric in an arbitrary gauge, 
then the gauge vectors £C(i)[7]j £^?(2)[7]7 etc., induce the 
transformation 

ft ^ ft' = ft + Aft 

= ft + e£e,,,.g+K(4(.,+4a,ff) 
-l-e'i;^,,,ft(i)-f ... (CI) 

If ft' is to satisfy the gauge condition [ft'] , then ^ must 
satisfy Lp[Aft] = — Lp[ft]. After a trivial calculation, this 
equation becomes 



^f-ne("„) = -si:ift«]-s^x"[ft(^)] 



n>0 



(C2) 

Assuming that this equation is solved for arbitrary 7, 
we can equate coefficients of powers of e, leading to a 
sequence of wave equations of the form 



(")' 



(C3) 



where is a functional of ^(1) , ^(n^i) and 

ft*^^', ft*^"^. I seek a solution in the region D, described 
in Sec. HIE. The formal solution reads 



on 



r\'Cin)~^?n^\'G%']dS'^ , (C4) 



Lorenz gauge in Fermi coordinates centered on 7: 



m 



1 i(i'i) 

2^a 



(C5) 



where is the mass dipole of the small body, 5* is its 
spin, and A'^''^'^ and Ca^'°'' arc irreducible pieces of ft'-^^ 
specifically 



c. 



(1.0) _ ,,(1,0) 



(C6) 
(C7) 



They are given in terms of tail integrals in Table I. 

The worldline is defined to be that of the body if 
vanishes for aU time. Given that A'P(t = 0) = dtM^{t = 
0) = 0, this is possible if and only if the right-hand side 
of Eq. (C5) vanishes. If, for simplicity, I neglect the Pa- 
papctrou spin term, then the first-order acceleration of 7 
must be 



a(i) - lini _ / !^ft(i) dfl - aft'^^ 



lim — 

r->-0 47r 

lim — 

j-->0 47r 



2'-'a'l-Ett 



i9ft('^ 

2^^"'Ett 



dth%l)dn 



dth 



(1) 



(C8) 
(C9) 
(CIO) 



where it is understood that explicit appearances of the 
acceleration are to be set to zero on the right-hand side. 
The first equality follows directly from Eq. (C5) and the 
definitions of ii^'^' and (7^^'°'. The second and third 
equalities follow from the STF decomposition of ft*-^' and 
the integral identities (D17)-(D19). The form of the force 
in the second line is the method of regularization used by 
Quinn and Wald [26] ; the form in the third line is used to 
derive a gauge-invariant equation of motion, as was was 
first noted by Gralla [70]. 

Now, suppose that we had not chosen a worldline for 
which the mass dipole vanishes, but instead had chosen 
some "nearby" worldline. Then Eq. (C5) provides the 
relationship between the acceleration of that worldline, 
the mass dipole relative to it, and the first-order metric 
perturbations (I again neglect spin for simplicity). The 

mass dipole is given by Mi = ^ limr_>o / r'^h^'^^riidQ, 
which has the covariant form 



where Gaa' is the retarded Green's function for the vec- 
tor wave-equation. From this we see that the Lorenz 
gauge condition can be adopted to any desired order of 
accuracy, given the existence of self-consistent data on 
the worldtube F of asymptotically small radius. I leave 
the question of that data's existence to future work. 

I now turn to the question of how the worldline trans- 
forms under a gauge transformation. I begin with the 
equation of motion derived in Ref. [28], written in the 



— lim 

87r r->0 



(Cll) 



where a primed index corresponds to a point on the 
worldline. Note that the parallel propagator does not 
interfere with the angle-averaging, because in Fermi co- 
ordinates, gp, = 5p + 0{e, r^). One can also rewrite the 
first-order-mctric-pcrturbation terms in Eq. (C5) using 
the form given in Eq. (CIO). We then have Eq. (C5) in 
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the covariant form 

-JTiaiV. (C12) 

Now consider a gauge transformation generated by 
eC'^^ [7] + jE^^*-^' [7] + where ^^^^ is bounded as r ^ 0, 
and Sp'^ diverges as 1/r. More specificaUy, I assume 
the expansions ^(i) = ^(^-O) (i, 6*^) + 0{r) and C^^) = 
+ 0(1).^° This transformation preserves 
the presumed form of the outer expansion, both in pow- 
ers of e and in powers of r. According to Eqs. (11)-(12), 
the metric perturbations transform as 

/^|,l^^/^i,l^+2e£.), (C13) 

Using the results for h'^^\ the effect of this transformation 
on IS given by 

h'^^^h^^-^-^n^e'+0[r-^). (C15) 

The order-l/r^ term arises from ft.|ii';p^^^'''' in the gauge 
transformation. On the right-hand side of Eq. (C12), the 
metric-perturbation terms transform as 

(C16) 

The only remaining term in the equation is mal^\ If we 
extend the acceleration off the worldline in any smooth 
manner, then it defines a vector field that transforms as 
a" — > a" +e£^(i,a" + .... Since a'-^^ = 0, this means that 
— it is invariant under a gauge transformation. 
From these results, we find that the left- and right- 
hand sides of Eq. (CI 2) transform in the same way: 

LHS/RHS ^ LHS/RHS 

(C17) 



The dependence on 7 appears in the form of dependence on 
proper time t. Each term could in addition depend on the accel- 
eration, but such dependence would not affect the result. 



Therefore, Eq. (C12) provides a gauge-invariant relation- 
ship between the acceleration of a chosen fixed worldline, 
the mass dipole of the body relative to that worldline, 
and the first-order metric perturbations. So suppose that 
we begin in the Lorenz gauge, and we choose the fixed 
worldline 7 such that the mass dipole vanishes relative 
to it. Then in some other gauge, the mass dipole will no 
longer vanish relative to 7, and we must adopt a differ- 
ent, nearby fixed worldline 7'. If the mass dipole is to 
vanish relative to 7', then the acceleration of that new 

worldline must be given by = eca + o{e), where 

(C18) 

Hence, this is a covariant and gauge-invariant form of the 
first-order acceleration. (By that I mean the equation is 
valid in any gauge, not that the value of the acceleration 
is the same in every gauge; under a gauge transformation, 
a new fixed worldline is adopted, and the value of the 
acceleration on the new worldline is related to that on 
the old worldline according to Eq. (C17).) An argument 
of this form was first presented by Gralla [70] for the case 
of a regular expansion of the worldline; it is now extended 
to the case of a fixed-worldline expansion. 

Appendix D: STF multipole decompositions 

This appendix briefly reviews the use of STF decom- 
positions and collects several useful formulas. Refer to 
Ref. [52, 71, 72] for thorough reviews. All formulas in 
this section are either taken directly from Refs. [71] and 
[72] or are easily derivable from formulas therein. 

Any Cartesian tensor field depending on two angles 9^ 
spanning a sphere can be expanded in a unique decom- 
position in terms of symmetric trace-free tensors. Such a 
decomposition is equivalent to a decomposition in terms 
of tensorial harmonics, but it is sometimes more conve- 
nient. It begins with the fact that the angular depen- 
dence of a Cartesian tensor Ts{9^) can be expanded in 
a series of the form 

Tsie^) ^Y.^s(L)n\ (Dl) 

where S and L denote multi-indices S = ii-.-is and 
L = ji...j£, angular brackets denote an STF combination 
of indices, n° is a Cartesian unit vector, = n^^ . . . n^', 
and = n''^^ . This is entirely equivalent to an expan- 
sion in spherical harmonics. Each coefficient Ts^^,^ can 
be found from the formula 

TsiD = J Ts{e^)nLdn, (D2) 

where !! is a double factorial (defined by x!! = x{x — 
2)...l). These coefficients can then be decomposed into 
irreducible STF tensors. For example, for s = 1, we have 

Ta(L) = f^at'^ + e\i.,fl°!^)^ + Sai^.f^ZW , (D3) 



30 











I 

~ 1 + 






_ 21' 




rpi-1) 






^ 2£- 


h 1 



, P<3 



where the T*^") 's are STF tensors given by 

(D4) 
(D5) 

(D6) 

Similarly, for a symmetric tensor Ts with s = 2, we have 

T^^L) = STF STF (e\,,fl;i\ + 5„,ri°li 

J- A 7^{^1) \ X X rp( — 2) 

T "(jjft Ut-l-'- pL~2 ^ L~2 



where 



{L+2), 

2£ 



(0) 



e + 2 L+i 

U{21~ 1 
(£+l)(2^ + 3) 



STF(T( 



jif) L-1 



(D7) 

(D8) 
(D9) 

(DIO) 



_ 2(£-l)(2£-l) 



-L-l 



T, 



(-2) ^ 



(^+i)(2£+TrL^!^^^<^-^>^'^-^"'-^''^' 

(Dll) 

2^-3, 



L-2 



2£+ 1 



(D12) 
(D13) 



These decompositions are equivalent to the formulas for 
addition of angular momenta, J = S + L, which results 
in terms with angular momentum £ — s < j <^ + s; the 
superscript labels (±n) in these formulas indicate by how 
much each term's angular momentum differs from (.. 

By substituting Eqs. (D3) and (D7) into Eq. (Dl), we 
find that a scalar, a Cartesian 3-vector, and the symmet- 
ric part of a rank-2 Cartesian 3-tensor can be decomposed 
as, respectively. 



T, 



(ab)\ 



e>o 
e>i 



(D14) 



(D15) 



e>o 



L-iri 



abL~2n 



L-2 I Ji f ^ L-2 

' If 



+ e.-' [aJ-b)jL-2ni 



l>2 



(D16) 



Each term in these decompositions is algebraically inde- 
pendent of all the other terms. 

The unit vector rii satisfies the following integral iden- 
tities: 



fiLdVL = if ^ > 0, 
riLdVl = if ^ is odd. 



(D17) 
(D18) 



n^dn = ^n ^^''Z''.^l\T' if ^ is even, (D19) 



(^ + 1) 



II 



where the curly braces indicate the smallest set of per- 
mutations of indices that make the result symmetric. For 

example, 5{abnc} = (Jab^c + <^fcc"a + Scam- 



Appendix E: Linear perturbation due to a point 
particle 

In this appendix, I present the solution to the wave 
equation with a point particle source. I also present the 
Detweiler- Whiting decomposition [59] of the solution into 
its singular and regular pieces. 



1. Solution to the wave equation 

The solution to the wave equation is^^ 

^eU - 2m /" G„;3a'/J'(2w"'u^' + (El) 

I seek an expansion of this equation in Fermi normal co- 
ordinates, in the case that x is near to a point on 7. The 
domain of integration can be split into two: the points 
in the convex normal neighbourhood M{x) — that is, the 
points that arc connected to x by a unique geodesic — 
and the points in the complement oi M{x). In the con- 
vex normal neighbourhood, the Green's function can be 
decomposed as 

Gapa'p- = Ual3a'P'S{a{x, x')) + Vafia' P'0{-a{x, x')). 

(E2) 

After performing a change of variables using dt' = — 
to evaluate the delta function, the metric perturbation 
becomes 

h^lp = —U^Po.'P'(2u"'uf'' + g"'^') + h'^l{u). (E3) 

^ict 

where primed indices now refer to the point on 7 con- 
nected to a; by a null geodesic; this point is given by 



This is not the standard form presented in, e.g., Ref. [54]; how- 
ever, one can easily show that it is equivalent to that form by 
using the Green's-function identities in Appendix A of Ref. [28]. 



31 



z°'{u), where u is the retarded tinie; rj-ot is the retarded 
distance between x and z°'{u)^ given by aa'U°' ; and the 
tail integral is given by 

Jt< 

r*^ , , , , 

+ 2m Gaf}a'fi'{2u°' +5" )dt' 

J —00 

= 2m/ Gc,/ja'0'(2M"'u'^' (E4) 



coordinates: 



where is the first time at which the worldline enters 
M{x). The two formulas for the tail are equivalent be- 
cause the upper limit of integration u~ = u — 0+ falls 
short of the past light cone of x, avoiding the divergent 
behavior of the Green's function there. 

The first term in h!^^\ sometimes called the "direct 
term" , can be expanded in powers of r using the follow- 
ing: the near-coincidence expansion [/q^q'^/ = (1-t- 

0(r^)); the transformation between riot and the Fermi 
radial distance r, given by 

T-cct = r[l + ^rai-n} — ^r^aiajn^^ — ir^daw" — ^r^hin} 



(E5) 



and the coordinate expansion of the parallel-propagators, 
obtained from the formula g'^ = 6^, where the re- 
tarded tetrad bf is given in terms of the Fermi tetrad in 
Eqs. (226)-(227) of Rcf. [54], and the coordinate expan- 
sion of the Fermi tetrad is given in Eqs. (123)-(124) of 
the same reference. The tail integral can be similarly ex- 
panded as follows: noting that u = t — r + 0(r^) , we can 
expand h^^'^u) about t as h^^^'^it) - rdth^^'\t) + ...] each 
term can then be expanded using the near-coincidence 
expansions V^^^ = gj^ g^^^ R"" ^" s" + 0{r) and 
h^fit) = g'^g^h'^l + rh'^%n') + 0(r2), where barred 
indices correspond to the point x = z{t), connected to 
a; by a spatial geodesic perpendicular to 7, and h^^^- is 
given by 



f^0, = 2™ / V^Ga^„,^,(2«"'u^' + g'^'^')dt'. (E6) 

'J —00 

This yields the expansion 

hlfiu) = g2g^^{hlf+rhlfy-imr£^-^)+0{r') . (E7) 

As with the direct part, the final coordinate expansion is 
found by substituting g^ = e"e^. 

Combining the expansions of the direct and tail parts 
of the perturbation, we arrive at the expansion in Fermi 



h')^l^ = —(1 -I- |rain* -f- ^r'^a^ajn''^ - ^r^daw" 
+ ^r^a^n' + ^r^Eiju'^) + (1 + 2rain')/i^^' 



2m, 



tail 



(E8) 



^Eta = 4maa - ^mrRoiajn''^ + 2mr£ain^ - 2mrda 

+ (1 + ra^n^Kf + rh^n^ + 0{r^) , (E9) 

+ ir^dj7i* - \r'^£ijn'-^)6ab + ^mraaab 

+ 0{r^). (ElO) 

As the final step, each of these terms is decomposed into 
irreducible STF pieces using the formulas (Dl), (D3), 
and (D7), to yield 

h%l^ = ^ + + 3m,a,ri* + r[Amaia' + if '^^n' 

+ m (|a(,a,) + f n'^] + 0{r''), (Ell) 
= C^o) + r(S(i'i)n. - 2md. + 

+ e„/i)f '^^n* + lmea,,Bih''') + 0{r^), (E12) 

- ^mSijh'^] + ^m£'i^^nh)i + 4nia(^aab) 

9 ''l-'^ab -r J^abi " ^ i'' b)j 



;^'%}+0(r2). 



(E13) 



where the uppercase script tensors are specified in Ta- 
ble I. The naming convention for those tensors follows 
that in Eqs. (D14)-(D16). 



2. Singular and regular pieces 

The Detweiler- Whiting singular field is given by 
/if^ = 2my' Gf^„,^,(2^."'u''' + g'''^')dt', (E14) 

where is the singular Green's function. Using the 
Hadamard decomposition 

Gf^a'^' = ^U^pCp'Sia) - iFa/Sa'/S'^M, (E15) 

we can write this as 

+ —Uaf!a"0"{2u""uP"+g''"^") 
^adv 

- 2m /" V^p^piu^uf^ + y^'^)dt (E16) 
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TABLE I: Symmetric trace-free tensors in the first-order met- 
ric perturbation in the buffer region, written in terms of the 
electric-type tidal field £ab, the acceleration ai, and the tail 
of the perturbation. 



^(1,0) 


I tail 
— "OO 




= /in"^'' -1- ma„ 
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^(1,0) 


L tail 
— ^{ab) 




J tail 1 n J tail , 2 
= /^OOa + ^ff'OO Cia + -^maa 




1 Ltail rzj 1 1 rtail z 


1) 

'^ab 




f)(l,l) 


— If i)c/ttail 1 Ltail \ 


£-(1,1) 


1 jcbc rtail 1 2 _ ■ 
= 3O "6ca + ^maa 


^(1,1) 


L tail 
— "(ai)c) 


^ a 


— 3 rijrtail 


f(i,i) 

^ab 





where primed indices now refer to the retarded point 
x' = z{u), where u is the retarded time; double-primed 
indices refer to the advanced point x" = z{v), where v 
is advanced time; radv is the advanced distance between 
X and 2:"(w), given by — o-q'/u" ; barred indices refer to 
points in the segment of the worldline between z{u) and 
z{v). The first term in Eq. (E16) can be read off from 
the calculation of the retarded field. The other terms are 
expanded using the identities v — u + 2r + O(r^) and 
radv = ?'rct(l + ^r^aiu'^); see Ref. [54] for details (though 
the expansion therein is for /if^.^, rather than /if^). The 
final result is 



Appendix F: Metric of a tidally perturbed black hole 

In this appendix, I present some general results for per- 
turbations of Schwarzschild in a light cone gauge. Over 
the course of the calculation, I highlight the restrictions 
that must be imposed on the metric in order to arrive 
at the usual result for a tidally perturbed black hole. In 
the final section of the appendix, I present that metric, 
along with its expansion in the buffer region. The no- 
tation and definitions in the first section mostly follows 
that of Ref. [641. 



1. Metric expansion, perturbation equations, and 
gauge condition 

The exact metric g is expanded as gi{X, e) = gsiX) + 
H{X,e), where H{X,e) = eH<^'^^X) + e^H'^^^X) + 
and X" = ([/, i?, Q^) are (scaled) retarded Eddington- 
Finkelstein coordinates adapted to the background met- 
ric gB, where R = R/e. As described in Sec. HIE, the 
terms in the inner expansion of g must satisfy the se- 
quence of equations (60)-(62), which I rewrite here: 



(1)1 - 



.(1) 



fiu r 



.9b J, 
(o)m!^[h-(i)1 



[9b] 



(Fl) 
(F2) 

)] 

(F3) 



2m 

+ f + 0(r2), 

2m „ 



{a'^b)i 



(E17) 
(E18) 



(E19) 



The regular field could be calculated from the regu- 
lar Green's function. But it is more straightforwardly 
calculated using h^^ = /i 



(1) 

a/3 



h^p. The result is 



/i,^ = i(i'°)+ri|^^^)n' + 0(r2). 



(7(1,0) +,,(5(1,1),, 

+ 0(r-^), 
6abk^^''^ 



(1,1) 



EJ-(l-l) • 

abi 



i f(l,l) i 



(E20) 

„■) 

(E21) 
(E22) 



In these equations, G^""* and S^C'f^^ consist of the terms 
in G/ and S'^Gi that contain n derivatives with respect 
to U. The first equation, (Fl), is the ordinary Einstein 
equation for gs, except that all derivatives with respect 
to U are removed. As the solution to this equation, I 
take the Schwarzschild metric 



gs = -f{U, R)dU^ - 2dUdR + R^dVt^, 



(F4) 



where / = 1 — '^'^^^^ , and M (U) is the Bondi mass of gB 
at time U divided by the initial mass. The dependence 
on U can not be determined at this stage, because time- 
derivatives appear only in the higher-order equations. 

Rather than fully solving the perturbation equations 
(F2) and (F3), I will solve only certain parts of them, 
in order to pinpoint several key points about the general 
solution. First, I adopt the light cone gauge. This gauge 



choice consists of setting H^^^ = H'^'^^ = H^^j^ = 0, which 
preserves the geometrical meaning of the retarded coor- 
dinates in the perturbed spacetime: U remains constant 
on each outgoing light cone, and R remains an affine pa- 
rameter on outgoing light rays. Second, as a boundary 
condition, I insist that the perturbations must be regular 
on the event horizon. 



(") 
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Because of the spherical symmetry of the background, 
it is convenient to expand the perturbations in tensorial 
harmonics: 



h7i 



^2 ^AB]^ 



tm 
AB 



(F5) 
(F6) 

(F7) 



where I have spHt the coordinates into the two sets = 
{U,R) and e^, the various harmonics will be defined 
below, and the coefficients of the harmonics are functions 
of U and R. In the context of this expansion, the light 



(n)£n 



cone gauge is imposed by setting 



P. 



RR 



RA 



0. 



I define the various harmonics as follows: The scalar 
functions are the usual orthonormal spherical 

harmonics, which satisfy [n^^DADg + i{i + = 
0, where ^ab is the metric of a unit 2-sphere, and Da 
is the covariant derivative compatible with ^ab- The 
even-parity vector harmonics Y^™ and odd-parity vector 
harmonics X^^ are defined as 



A 



DaY 



A 



-ea^DbY 



(F8) 



where cab is the Levi-Civita tensor on the unit two- 
sphere. The even-parity and odd-parity tensor harmonics 
Y^'g and X^ are defined as 



Yi^ = [DaDb + + mAB\ F'™, (F9) 
^Tb = -1 {^a'^Db + eB^DA) DcY'^. (FIO) 



In this appendix, I will forgo any discussion of the odd- 
parity terms, since the even-parity terms are sufficient 
for my purpose of delineating the types of restrictions 
required to arrive at the usual form of a tidally perturbed 
metric. 

In order to determine the effect of a gauge trans- 

^(1) _ 



formation, I write an even-parity gauge vector 

■^A 



^(1) 



^(1) 
^A 



m\rzm 



'Y, 



(Fll) 



This vector has the following first-order effects: 
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(F12) 
(FIB) 
(F14) 
(FI5) 
(F16) 
(F17) 
(FI8) 



where for simplicity I have omitted the harmonic labels 
£m. I neglect [/-derivatives in the transformation, since 
they are second-order effects in the present scheme. 

Now, due to the spherical symmetry of the background 
metric, each mode in the harmonic expansion decouples 
from all the others in 5Gi. Likewise, the even- and odd- 
parity sectors decouple. I write the even-parity terms in 
SGf'> as 

Qtm = J SGf^^Y^'^dVl, (F19) 

j SGf^^^nABY'"'dn, (F21) 

Explicitly, these gauge-invariant quantities are given 
by [64] 



92 2 d ~ f d ~ 2 d ~ 
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(F23) 
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- ^<i±il±ift„ + Xp™ + JLk. (F24) 

{R-M)f d ~ f d ~ uR + 4M~ 

U = = —K — — -riJU H = t^UU 
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(F25) 



(F26) 
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(F27) 



2{R-M) g {R-M)f dj^^^ 
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i?2 dR i?2 

+ 1)^2 _ 2^lMR - 41/2 ^ 



Prr, (F28) 
(F29) 



where /i = 1) — 2, and I have introduced the gauge- 
invariant combinations 

2M 2Mf d 
Puu = Puu - -^Ju + -^Jr - Mf-~-G, (F30) 
R R oR 

d ,~ , 2M d , , 

Pur = Pur - T^iRJu) + -^Jr - M-^G, (F31) 
uR It oR 

d ~ - (92 ~ d 

= Prr - 2—^(RJr) + i?2 — _G + 2R—^G, 
dR^ ' 9i?2 

(F32) 



K = K + 2Ju- 2fjR + Rf-^G + Ml + 1)G. 

dR 



(F33) 



I have omitted the indices {n)£m for simphcity. 



First-order solution 



The first-order equation reads 6G^P°'^[H^^''] — 

—G^'plgs]- The source term in this equation has a single 
non- vanishing component, 



g[''^^[9b] 



2 dM 



(F34) 



Thus, in the notation introduced above, the RR, H ^ Q 
equation reads Q^^qq = '^V^P~^^if^ while all the other 
equations are source- free. For £ > 2, the equations can be 
solved for arbitrary i. However, because various quan- 
tities are defined only for £ > 2, the equations for the 
low multipoles £ = and £ = 1 must be dealt with in- 
dividually. I will write undetermined functions of U as 

4")""(f/). 

I begin by solving the £ = equations. For £ ^ 0, 
the quantities Ja, G, Q^, and are undefined. So 
the only equations are Qfo^ooi^''^^] = Qfofooi^'-^^] 
QmooiH^'^] = and QgUHi^^] = ^V^R-'j§, in 



which Ja and G are set to zero. Given the gauge condi- 
tion, the only functions appearing in these equations are 
i^(i)oo and P^^j}°°. I first solve Q([,)^oq[-^^^^] = 0' ^^^ich 



reads explicitly 



dR^ 



R OR 



= 0. The solution 



to this equation is 



^(1)00^^(1)00^4^(1)00^ 

^ R ^ 



(F35) 



Substituting this into Q^J^qqIH''^^] = yields an equation 

d(1)0 

uu 



for P^l}™ that can be readily solved to find 
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(F36) 



Substituting these results into (5(Q-|QQ[i/^^''] and 
QfofooiH^^^], we find that both quantities are 
identically zero. Hence, from the equation 

Qfo^olH'^^^] = 4:^R-^^ we can conclude 



dM 

lu 



0; 



(F37) 



that is, the mass of the internal background is constant. 
The functions A'j}^^^ , k = 1,2,3, can be be determined 
only by solving the second-order EFE. 

Next, I proceed to the £ = 1 equations. For £ — 
1, the quantities G and are undefined, and the 
field equations read QI^^.JH^'^] = Qlo)iJH^'^] = 
'3(o)im[-^'^^] ~ which G is set to zero. Solving 

Qyonm = yields 
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^ R 



(F38) 



Substituting this into QYn)^,-^ = and solving then yields 
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Substituting these results into = and solving 

then yields 
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(F40) 

Two of the remaining equations, Q^^^^ = = Q^fim 
fixes several of the free functions in these solutions: 



a: 
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and A 



(l)l7n 



0. The final equation. 



^(o)im ~ yields no further information. Putting these 
results together, we find 



3(1) Ir^ 
UU 



(l)l^ 
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l)lm 



{l)lm 



P2 " R2 
J_^(l)lm ^ J_^(l)l'i 
2P ^ P2 3 
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(F41) 
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Finally, I proceed to the £ > 2 equations. As with 

)UU 

'{o)en 



0, 1, the equation = can be immediately 



solved to find 



^ R ^ 



(F44) 
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Using this result and Qfrs 



{0)lm 



0; wc can express Fjjjj 



in terms of G'^^-'^'" and J^''^"^; next, we can use Q'^q^^ 



to express G^^-'^™ in terms of J^''^" 
ing these results into Q^q^^,; 
to find 



finally, substitut- 



0, we can solve for J, 



(l)«m 
U 



A, 



l)^m 



2M 



(F45) 

where 2-F1 is a hypergeometric function. The term next 
■ diverges at the unperturbed event horizon, R = 
2Af, violating the boundary condition. Therefore, we 
have Ag"'^'^™ = 0.^^ Next, I make my first restriction 
of the solution; if J^^-*^™ is expressed in terms of i? = 
eR, then the term next to ^W^™ ^^gj^^ves as e~^. Since 
H^^^ is accompanied by a factor of e in the metric, this 
term behaves as e^^^. If I assume that R ^ r, where 
r is, for example, the Fermi radial coordinate centered 
on the body's worldline in the external spacetime, then 
these negative powers of e would also appear in the outer 
expansion. By assumption, no negative powers do appear 



in the outer expansion; therefore, I set 



(l)£m 



0. So 



Recall that we had expressed P^^"' and G'-^'^"' in 
terms of J^'^^^^ . With j'jj^^^^ determined, we now have 
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Substituting these results into Q^'q^^^ = and then 
Qoi^im = we determine 



'(0) 



(i)£ti 



2A 



(l)£n 



1)' 



(l)£m 



(F49) 



The sole remaining equation, Q(^i^£^ — 0, yields no new 
information. Hence, the first-order calculation is now 



complete. We have found that = 0; the £ = modes 



in are given by Eqs. (F35) and (F36); the £ = 1 

modes are given by Eqs. (F41)~(F43); and the £ > 2 
modes are given by 



P, 



(l)fn 



UU 



f^{Ai''>'--2£i£+l)Ai'^'-), (F50) 



J, 



(l)fon 



M 
R^ 

2i? ' R^ ^ ^ ' 



1 



R 



Q(l)ln 



-^A 
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2A 



{l)im 



£{£ + !) 
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These results can be simplified by a refinement of the 
lightcone gauge. For £ = 0, the function A'^^'^^ can be 
removed via 

.(i)oo^rn (F54) 



c(i)oo 



f(l)00 _ .<;{1)00 _ 1 ,(1)00 
?(7 ^ -I^R 2^2 ' 



leaving 
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A 



(1)00 



-.A 
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1(1)00 



(F55) 
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Although this does not exhaust the residual freedom 
within the lightcone gauge, since ^r{U) is arbitrary, 
the remaining freedom cannot be used to remove either 



(1)00 



A 



(1)00 



However, if we were to transform out 



of the lightcone gauge, A^"'^-'™ could be removed, leav- 
ing only a 1/R mass monopole term, corresponding to a 
time-dependent correction to the mass. 



For £ > 1, we can removed all A 



(l)fn 



Via 



f(i)i™ _ 14(1)1'- 

— ~2"^2 



1 .(l)lm 

R ^ 



^(l)lm 
f(l)lm 



2M 



-A 



(l)lm 



(l)l?ri 



and 
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1 4(1)^™ 
'2^2 



2M 



1 4(l)lm 
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{l)em 



{l)lm 
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(F58) 
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(F61) 
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The term next to Ag^'^™ corresponds to an induced multipolc 
moment. The fact that it must vanish agrees with the no-hair 
theorem. 



This exhausts the residual freedom in the gauge condi- 
tion. 

In order to arrive at the usual form of a tidally per- 
turbed metric, we must go beyond these gauge refine- 
ments by setting the entirety of H^^^ to zero. That is. 
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we must choose A^"'^^"" = ^3^''°° = 0. If the odd-parity 
calculation had been performed, we would find that £ = 1 
terms, corresponding to time-dependent spin terms, must 
be set to zero. (Recall also that I have restricted the pos- 
sible perturbations by disallowing terms that would con- 
tain negative powers of e in the unsealed coordinates.) 
Hence, we can conclude that to arrive at the usual met- 
ric of a tidally perturbed black hole, we must restrict 
the perturbation by setting numerous functions to zero, 
without any evident justification. 



3. Second-order solution 

With set to zero, and with M determined to be 
a constant, the second-order EFE becomes the homoge- 
neous, linear equation 



(F64) 



This equation is solved in the same manner as the first- 
order equation. The calculation of the low multipoles 
£ = and £ ^ 1 proceeds just as at first order, yielding 



and so 

j(2)2m _ 1 .(2)2m 4M-I-4R 

A (2) 2m 

^ 2R ^ m 

- A^^^^'^'R^f. (F72) 

As at first order, this determines P^]^™' and G^^-'^™: 
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uu 



^2^2 4 



(2)2r^ 
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(2)2m 



3^(7?-4M)yl^^^^", (F73) 
4 



Q(2)im ^ _lj^mm ^ ^242)2™ ^ ^(2)2™_ ^^^^^ 

R 



As at first order, we determine Ag^-*^™ and Ay^''^™ from 
QYoUm = and then Q/^{^. = 0, which yield 



(0)£n 

(2)2)71 _ _ 1 ^{2)2m _ 27Vf2^(2)2m 



3^1 
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UArA 



,r2 x(2)2m 
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(F76) 



Putting these results together, we have the solution 



P, 



_ .(2)00 _ M (2)00 1 



(2)00 



(2)00 



_^(2)00 ^ ^(2)00 



1 



(2)00 



and 
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UU 



J2)lr 
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For t > 2, the calculation proceeds just as at first 
order, up until the point marked by Eq. (F45). When 
the term that diverges on the event horizon is removed, 
the analogue of that equation reads 



J, 



{2)tm _ 1 <(2)toi 4M + y,R .{2)ln 

2R ^ R^ 



2M\ 



~R!) 
(F70) 



At first order, this solution was simplified by setting the 
coefficient of R^ to zero, because it led to negative powers 
of e when written in terms of R. However, at second 
order, this term will be multiplied by in the metric, 
so it will scale as e^~^. Hence, the term is acceptable for 
£ ~2^ but must be set to zero for £ > 2. I will hence deal 
with these two cases separately. 
For £ = 2, we have 



2FA 2 



Rf 



(F71) 
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Ju --^^2 +^11 + ^1 A, 
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G(2)2™ = _iA(2)2'" +R^(l 
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(F79) 

(F80) 



For £ > 2, the calculation is proceeds just as at first 
order, leading to the solution 



= -^A'" + ^^(^ + (FBI) 

R'^ R'^ 
= --Ia^^)^'" + ^^ + ^^^ Af)^"\ (F82) 



Q{2)lr, 



£{£ + !) 



(2)hn 



-^A 

R 



{2)e.m 
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Just as at first order, with an appropriate gauge re- 
finement, we can remove the functions Aj^''^", A[,^'^™, 
^(2)im^ j^(2)irn^ ^^^^ ^(2)im ^ ^^^^^ ^ ^ 2, thereby ex- 

hausting the freedom within the lightcone gauge. In or- 
der to arrive at the usual form of the tidally perturbed 
black hole metric, we must then set Ap-*"*^ = Ag'^-'"*^ = 0. 
(And I again remind the reader that I have disallowed 
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terms that would possess a negative power of e when ex- 
pressed in unsealed coordinates.) This leaves only one 



undetermined function: A 



{2)lm 



Although I do not show 



the odd-parity calculation, it yields an analogous result: 
to yield the usual form of the metric, all but one of 
the possible undetermined functions must be set to zero, 
leaving a single quadrupole term. After imposing these 
restrictions, the only non-vanishing components of the 
metric perturbation arc 

^3i?2/'4'^'"5^"", (F85) 
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(F86) 
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4. Tidally perturbed black hole metric and its 
expansion in the buffer region 

After the metric perturbations have been restricted as 
described in the preceding two sections, the inner expan- 
sion has the following form: 

^f + e^Y. 3i?V'4'^""i""" + (F88) 
-R e^R^fJ2A^^^"'YA"' 

m 



giuu 



giuA 



giAB = R 



AB 



e'R^ ( 1 - 
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;2)2m^2' 



Y, 



AB 
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along with the exact results giuR. = ~1 and giRR ~ 
giRA = 0. Although it is written in a slightly different 
form, this is the usual metric of a tidally perturbed black 
hole. It is characterized by (i) having only quadrupole 
perturbations, and (ii) those perturbations scaling as R^ 
for large R. 

To write this metric in terms of a pair of tidal fields 
Sab and Bab, I follow Appendix A of Ref. [73] ; I also make 
use of notation and simple identities from the Sec. 3.3.7 
of Ref. [54]. First, I define the quantities 



5 ^ ' 



m 



(F91) 
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Next, I define the derived quantities 
^A = \DaS* , 



and 



£\b = {DADB + 3nAB)£*, 

B*AB = -5 {^a^'Db + eB^DA) DcB* 



(F93) 
(F94) 

(F95) 
(F96) 



Using the definitions of f*, B\ yj™, Y^, X^/", and 
-^AB' '^^^ express the derived quantities as 



and 



£*ae 
B*A-= 

^AB 
B*AB 



E^. 



(2)2m-^2m 



|^^(2)2m^2m^ 
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m 

_3^5(2)2m^2™^ 



(F97) 
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(F99) 
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In terms of these quantities, we can write the metric as 
giuu ^-f ~e''R\f£* +0{e^), (FlOl) 



giuA = R 
giAB = 



^e^R'f (£* +B*) +0{e'') 



Uab \e'ie f 1 - ^ 1 f * 



AB 



\e'R^B*AB + Oie^) 



(F102) 



(F103) 



This is the usual form of the metric of a tidally perturbed 
black hole. 

Now, because £* and S* are, respectively, even- and 
odd-parity quadrupole terms, they can be written in 
terms of an STF decomposition: 

£* - fabfi^"") , B: = EacdBtn'^"''^ . (F104) 

By applying the definitions (F93) and (F95), we find 
£*A = n'XEabn", B*A = n'XeacdBtn^"-^ , (F105) 

and 

SXb = '^^A^B^ab + nABSab^'-''''^ , (F106) 

B*AB - n'xn'BeacdB^n' + n'^n^etcdBin^. (fio7) 

When written in more explicit form, these expressions 
agree with those in Ref. [54]. In order to arrive at these 
results, one requires the identities 



iAB = eabc^A^B^'^ 



^ Bob 
€A "r 



DADBn" = -n^^AB- 



(F108) 
(F109) 
(FllO) 
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Finally, I convert to Cartesian coordinates ([/, F"), 
adapting the identities of Sec. 3.3.7 of Ref. [54]. The 
result is 



(Fill) 

a, + ^e^R'fiS: + Bl) + 0{e'), (F112) 



giuu = -f - + 0(£^), 

giUa — --a I 3 



giab — Sab — ^ab — he^R^ 1 — 
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''ah 



Se^R^Bl, + 0{e^). 



(F113) 



where = £\nt Bl = B\nt £1, = £1b^^^^ , and 
i5*^ = B\gVt^Vt^ . When this is expanded in the buffer 
region, by rewriting it in terms of the unsealed radial 
function R and then re-expanding for small er, it becomes 



giuu 



giua 



2 A/ 

-1 + e— - R^£* + 4eMR£* 
R 

-0{e\eR^,R^) , 

+ lR\£: + Bl) ~ \MR{£1 + Bl) 

2 n3\ 



0{e,eR,R 



(F114) 

'*a) 

(F115) 



giab = Sab - Nab ^ ^^^(^ab + Kb) 

+ o(£^ei^^i^^) . 



(F116) 



In order to agree with the external Lorenz gauge, I 
switch to harmonic coordinates via the transformation 
Y" = X""' + eMN''', where TV"' = X"' /R' (note that. 



in the buffer region at first order in e, this transfor- 
mation to harmonic coordinates cannot be distinguished 
from a transformation to isotropic coordinates). I then 
switch from retarded coordinates to Fermi-type coordi- 
nates {T,X°-) via the transformation 

U = T-R~ 2eM{\nR+ \R^£^jN''^) 



(F117) 

X"' ^X" - iR^R^bocN'"' + IR^£SN'', (F118) 

where iV° = X°- /R. After performing these transforma- 
tions and decomposing the result into irreducible STF 
pieces, we arrive at 

2M 



giTT = -1 + e— + 'pMR£,jN^^ - R^£,jN'^ 

-^0(e^£i?^i^3) , 
giTa = 2eMR£a.,N' + ^eAIRea^MN''' 

2 n3^ 



(F119) 
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3"" ^aiK^j 
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giab = Sab{l + e— - leMR£,jN'^ ~ lR^£^,N''' 
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\eMR£,(^aNb)' - ^eMREab - lR^£ab 



leMR£,,Nab'' + lR^£,(aNl^ 



\eMRe^k(aBtN'' + 0{e\ eR\ R^) . (F121) 
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